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WAVE PROPAGATION ON EUCLIDEAN SURFACES WITH 
CONICAL SINGULARITIES. I: GEOMETRIC DIFFRACTION. 

G. AUSTIN FORD, ANDREW HASSELL, AND LUC HILLAIRET 


Abstract. We investigate the singularities of the trace of the half-wave group, 
Tre“*^^^, on Euclidean surfaces with conical singularities {X^g). We com¬ 
pute the leading-order singularity associated to periodic orbits with successive 
degenerate diffractions. This result extends the previous work of the third 
author [Hil05] and the two-dimensional case of the work of the first author 
and Wunsch [FW] as well as the seminal result of Duistermaat and Guillemin 
[DG75] in the smooth setting. As an intermediate step, we identify the wave 
propagators on X as singular Fourier integral operators associated to intersect¬ 
ing Lagrangian submanifolds, originally developed by Melrose and Uhlmann 
[MU79]. 
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0. Introduction 

In this article, we investigate the spectral geometry of Euclidean surfaces with 
conical singularities {X,g). We determine the precise microlocal structure of the 
half-wave propagator, near a ray that undergoes one or two degenerate 

diffractions. Using this, we compute the leading-order singularity of the trace of 
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the half-wave group, associated to an isolated periodic orbit undergoing 

two degenerate diffractions through cone points. For example, if the periodic orbit 
has length L and undergoes degenerate diffractions through two cone points at a 
distance b apart, we show that the associated wave trace singularity is 

( 0 . 1 ) 

0.1. Background. Spectral geometry typically aims at understanding the rela¬ 
tions between the spectrum of the Laplace operator on a Riemannian manifold and 
the geometry of the associated geodesic flow. These relations may be revealed by 
the study of wave propagation. For instance, the Poisson relation states that the 
trace of the wave propagator is smooth except possibly at the lengths of periodic 
orbits. Moreover, in a generic and smooth situation, the singularity that is brought 
to the wave trace by a particular periodic orbit can be fully understood and leads to 
the definition of the so-called wave-invariants (see [DG75]). These wave-invariants 
may then be used for instance in inverse spectral problems. They also serve as 
a particular motivation to study wave propagation on different kind of singular 
surfaces. We will focus on Euclidean surfaces with conical singularities since this 
general setting includes polygonal billiards and translation surfaces, both of which 
are very interesting and natural. 

The basic new feature of wave propagation on singular manifolds is the di¬ 
chotomy between waves that hit the singularity—that are then diffracted in all 
possible directions—and waves that miss the singularity and propagate according 
to the usual laws for smooth manifolds. This fact leads to the definition of the 
so-called geometric (or direct) front that consists of rays that miss the vertex and 
the diffracted front that consists of rays that hit the vertex and are reemitted in all 
possible directions. On a two dimensional cone, these two fronts share two rays in 
common that correspond to the limit of rays that nearly miss the cone point from 
one side or the other. In the literature, these two rays are called either “geomet¬ 
rically diffractive” [MW04] or “singular diffractive” [Hil05]. We will use here the 
former terminology. On a compact surface with conical singularities the situation 
becomes quickly complicated for a diffractive ray may hit successive conical points 
and experience new diffractions that may be singular and so on. These diffractive 
phenomena are established in the quite abundant literature on wave propagation on 
singular manifolds starting with Sommerfeld’s result for Euclidean sectors or cones 
[Som96]. Among the important milestones of this story are the studies by Cheeger 
and Taylor for cones of exact product-type [CT82a, CT82b] and by Melrose and 
Wunsch in the general case [MW04]. 

Over the years, there has been investigation of the impact of diffraction on the 
wave-trace. Eor instance, Wunsch showed in [Wun02] that singularities may appear 
at length of periodic diffractive orbits. For some periodic diffractive geodesics, the 
leading singularity is then computed in [Hil05] in the Euclidean case and in [EW] in 
a more general case (see also [BPSOO] for related results from a physics perspective). 
Both these results are built upon a precise description of the wave propagator that 
is microlocalized in the vicinity of given periodic (possibly diffractive) geodesic. 
However, none of these studies attempted to determine the precise microlocal na¬ 
ture of the propagator near a geometrically diffractive ray: in [EW], it is assumed 
that no geometric diffraction occurs (with a non-focusing assumption that would 
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be automatically satisfied in our case), while in [Hil05], it is assumed that the pe¬ 
riodic geodesic has at most one geometric diffraction. The main purpose of the 
present paper is to fill this gap, i.e., to give a precise microlocal description of the 
wave propagator near the geometric diffractive rays, on an ESCS. More precisely, 
we will identify the microlocalized propagator near a ray that undergoes one or 
two geometric diffractions as an element of the Melrose-Uhlmann class of singular 
Fourier Integral Operators ([MU79]), associated to either two, or four, Lagrangian 
submanifolds. One advantage of this identification is the ease of computing wave 
trace singularities, such as (0.1), using standard methods such as stationary phase. 

This is the first article in a planned series of three. In the second paper, we 
will show how to compute wave traces for any closed orbit on an ESCS (with any 
number of geometric diffractions). In the third paper, we will apply our results to 
inverse spectral results, specifically isospectral compactness in the class of ESCSs. 
To keep the length of the present paper within reasonable bounds, we restrict our 
attention here to at most two geometric diffractions. 

0.2. Cones and ESCSs. The Euclidean cone of cone angle a > 0 is the product 

def 

manifold Ca = (R+)r ^ equipped with the exact warped product metric 

The vertex of the cone p is the point where all (0, 6) are identified, and we will 

def 

denote by C° = \ {p} the cone without its vertex. Let us recall that the 

Euclidean distance on Ca between two points qi = (ri,9i) and <72 = (t 2,02) in 
polar coordinates is: 

{ dist(p,qi) = ri, 

dist(qi,q2) = Ti-br2, |6»2-6'2|>7r 

dist(qi,q' 2 ) = -v/d -b - 2 rir 2 cos( 6»2 - 6 >i), | 6»2 - 6 »i| < tt. 

A Euclidean surface X with conical singularities (denoted by ESCS in the sequel) 
is a singular Riemannian surface such that any point has a neighbourhood that is 
isometric either to a Euclidean ball in or to a ball centered at the vertex of some 
Euclidean cone Ca- 

Example 0.1. From any polygonal domain If in the plane we may generate an 
ESCS by taking two copies of the polygon, reflecting one of these copies across the 
y-axis, and identifying the corresponding sides. Starting from a square, we build in 
this way a surface that is topologically a sphere that is flat with four singularities 
of angle tt. 

Example 0.2. More generally, a surface that is obtained by gluing Euclidean 
polygons along their sides also is Euclidean with conical singularities. The surface 
of a cube is a ESCS that is topologically a sphere with 8 singularities of angle 

Let X be a Euclidean surface with conical singularities, and let P be the set of 

def 

its conical points. Define X° = X \ P. Let m be a smooth function that vanishes 
near the conical points. Using the Euclidean metric, one defines the gradient of u, 
Xu, and the action of the Laplacian on u, Am, as usual. The Laplace operator thus 
defined is not essentially self-adjoint. Among the possible self-adjoint extensions. 
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the most natural one is the Friedrichs extension that is associated with the Dirichlet 
energy quadratic form 

Q{u)= f \Vu\^dS, uGC^iX°), 

Jx 

where dS is the Euclidean area measure. Throughout the paper A will always 
define the Friedrichs extension of the Euclidean Laplace operator. By choice it is a 
non-negative operator. 

Writing □ = — A with Dt = j dt, the associated wave operator is then defined 

as 

x) = 0 

u( 0 ,x) = uo(x) 
dtu( 0 ,x) = uo{x) 


We will always take < ^ 0. The wave propagators that are associated with this 
wave equation may be defined through functional calculus and we denote them by: 


(0.4) 


W(t) 


i(tyA) 

7 ^ 


and W (t) cos (tyfK 


We will also use the half-wave propagator U{t) = exp{—itVX). 

Since singularities of solutions to the wave equation propagate with finite speed, 
the propagator W(t) can be understood by patching together local propagators 
that are defined either on the plane or on C^- As a first step it is therefore crucial 
to understand wave propagation on the flat cone Ca- 


0.3. The wave kernel on cones. It turns out that the wave kernel on Ca is ex¬ 
plicitly known (see [Som96,CT82a, CT82b,Fri81] for different ways of constructing 
this kernel — we describe these briefly at the beginning of Sections 2 and 3). Prop¬ 
agation of singularities for the wave equation on Ca is then described as follows. 
Using polar coordinates, we define on (0,oo) x T*C° x T*C° two Lagrangian sub¬ 
manifolds and A*^. For a > tt, these can be defined as follows. The geometric 
(or “main”) Lagrangian is 

(0.5a) = rl +r 2 — 2 rir 2 cos(di — 62 ) and \9i — ^ 2 ! ^ tt} , 

the diffractive Lagrangian is 

(0.5b) =7*{t2 = (ri+r 2 )"}, 

and their intersection is the singular set 

(0.5c) S = A°n{| 6 »i- 6 i 2 | = 7 r}. 

In the case a < tt, we choose an integer N such that aN > tt. Then we consider the 
A-fold covering map from to (7° induced by the natural map M/AaZ —>• K/aZ. 
As this is a local isometry, this induces a covering map T*C'^^ —)■ T*Ca- We define 
A^ to be the image of A^^ under this covering map. 

The terminology indicates that A^ corresponds to geometric, or non-diffractive 
geodesics (i.e., geodesics on Ca that avoid p) which carry the main singularity 
whereas A'^ corresponds to diffractive geodesics (i.e., concatenation of two rays 
emanating from p.) The singular set thus corresponds to diffractive geodesics that 
are limits of non-diffractive ones. We will refer to these as geometrically diffractive. 
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We will denote by the Lagrangian submanifolds obtained by restricting 
to =Fr > 0 where r is the dual variable to t. 

The explicit expression of the propagator implies, first, that singularities prop¬ 
agate according to A*^ U A°, and second, that away from the intersection E the 
propagator is a classical Fourier Integral Operator (FIO). Away from the inter¬ 
section E, the kernel of the half-wave propagator is given by the so-called 

Geometric Theory of Diffraction (see Appendix B). 

0.4. Main results. Our first result is a precise description of the kernel of the 
wave propagator on the cone Ca near the singular set E. It is actually a bit simpler 
to describe the result for the half-wave propagator whose Schwartz kernel 

we denote by Ua- 

We observe that (t*, <7i, 52 ) is in the projection of E on (0, 00 ) x C° x C° if and 
only if, in polar coordinates, we have r -|- = t* and 91 — 62 = ctt, e = ±1. Let 7 

be the parametrization by arclength of the geometrically diffractive geodesic that 
joins ql to normalized in such a way that 7 (—r^) = < 72 , 7 ( 0 ) = p, 7 (^ 1 ) = qi- 
Since the cone is flat, 7 can be extended to a local isometry that is defined on 

\ {(0,ej/), y > 0, e = ±1}. Using we can thus parametrize a neighbourhood 
of (QiiQI) ill C'° X C° by the product of two Euclidean balls in the first one 
centered at (rj',0) and the second one at (—?’|,0) (in Euclidean coordinates). 


Theorem 0.3. Let ql and ql he the extremities of a geometrically diffractive ge¬ 
odesic of length t* and diffraction angle eir (e = ±1). Locally near {t*,ql,ql) in 
(0, 00 ) X (7° X (7° the kernel Ua can be written as the following oscillatory integral: 

(0.6) Ua{t,qi,q2) = i 2 TT)~'^ f f qi, 92, s, w) dwds 

J s>0 J u;>0 

where (using for parametrization — i.e., gffxi,yi) = qi, geix 2 ,y 2 ) = Q 2 J 
(1) the phase (/>e is defined by 


((e{t,qi,q2,s,uj) = UJ 


\Jxl + {yi + seff 4- y/xff+Jy^^^s^ - t 


( 2 ) the amplitude aa,e is a classical symbol that is smooth in (t, gi, < 72 , s) and 
of order 1 in uj so that we have 


^ ^ —/i:(t, Ql, ^2; 


UJ 


1 -k 


k >0 


(3) In polar coordinates, we have at leading order 

( n ^ o ■ Sa{ 0 i- 92 ) 

aa,€{qi, 92 , S = 0,w) = -2me ■ 


sin 01 -I- sin 02 


(rir2)2 

where Sa is the (absolute) scattering matrix for the cone Ca- An explicit 
expression for Sa is given by (B.ll). 


From this expression we deduce the following corollary. 


Theorem 0.4. The half-wave propagator Ua{t) on the Euclidean cone Ca is in 
the Melrose-Uhlmann class/'"(A°, A^) of singular Fourier Integral operators. The 
order m is equal to 0 if t is regarded as a parameter, or —1/4 if t is regarded as a 
‘spatial’ variable. Similarly, the sine propagator'W(t) on the Euclidean cone Ca is 
in the Melrose-Uhlmann class /'"“^(A'^, A^) of singular Fourier Integral operators. 
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It can be noted that elements of this class are standard FIOs away from the 
intersection S so that this theorem doesn’t say anything new away from E. On 
the other hand, although the explicit expression of the propagator was already 
known near E, the fact that it belonged to the Melrose-Uhlmann class was not. It 
is also worth remarking that it may be possible to obtain the latter theorem by 
some brute computations starting from the explicit expression of the propagator. 
We propose a different method, the ‘moving conical point’ method, that exploits 
geometric features of wave propagation on cones. It has the advantage that the 
parameter s in (0.6) then has geometric significance: it is the distance by which 
the conical point is shifted. 

Remark 0.5. It is actually convenient to use the Riemannian metric to identify 
functions and half-densities. This amounts to multiply the oscillatory integral rep¬ 
resentation by the half-density \dq'dq \2 or \dtdq'dq\ 2 . 

Remark 0.6. Recall (or see Section 1) that in the Melrose-Uhlmann calculus, the 
order of the distribution on the first Lagrangian is ^-order less than on the 
second, A*^. This allows to recover the fact that the diffracted wave is ^-order 
smoother (in a Sobolev sense) than the direct wave (in two dimensions). 

The oscillatory integral representation of the preceding theorem has several in¬ 
teresting applications mainly because it allows one to compute simply the wave 
propagator on an ESCS when microlocalized near a geodesic with several geometric 
diffractions. We will illustrate this by obtaining, for a geodesic with two geometric 
diffractions in a row an oscillatory integral representation that fits into the class 
of singular FIO that is constructed in [MU79, Sections 7-10] and associated with 
a system of four intersecting Lagrangians. More precisely, consider a geodesic of 
length t between q and q' with two geometric diffractions at pi and P 2 ■ There are 
four types of nearby geodesics: 

(1) non-diffractive geodesics; 

(2) geodesics that are diffractive at pi but not at P 2 ; 

(3) non-diffractive geodesics at pi that diffract at p 2 ', and 

(4) geodesics that diffracts at pi and P 2 ■ 

Each type corresponds to a particular Lagrangian and these four Lagrangians form 
a intersecting system in the sense of [MU79]. 

Using the preceding theorem and standard stationary phase arguments we obtain 
the theorem. 

Theorem 0.7. Microlocally near a geodesic with two geometric diffractions, the 
half-wave propagator on a ESCS is in the Melrose-Uhlmann class of operators as¬ 
sociated with a system of four intersecting Lagrangians. 

We actually get much more accurate information since we can derive the principal 
symbol of the half-wave propagator on the twice diffracted front — see (4.15) and 
(4.16) . 

Finally we will use our new expression for C/q, to compute the contribution to 
the wave-trace of an isolated periodic geodesic with two geometric diffractions. 

Proposition 0.8. On a ESCS, the leading contribution to the wave trace of an 
isolated periodic diffractive orbit with two geometric diffractions is 

1 

4z7r^ 


• \/b{L — b) ■ {t — L — iO) ^. 
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This is perhaps the simplest setting for which neither [FW] nor [Hil05] applies. 
This proposition shows that such a geodesic creates in the wave-trace a singularity 
that is comparable to the singularity that is created in a smooth setting by an 
isolated periodic orbit. The singularity is ^ stronger than a diffractive geodesic 
with one non-geometric diffraction and | weaker than a cylinder of periodic orbits. 

With our new representation of the wave kernel, it should actually be possible to 
compute the full asymptotic expansion of the contribution to the wave-trace of any 
kind of periodic diffractive geodesic. This is a far-reaching generalization of results 
in [Hil05] and it leads to the possible computation of many wave-invariants. This 
opens new questions concerning inverse spectral problems in this kind of geometric 
setting which, we recall, includes Euclidean polygons. For instance it can be asked 
whether the full asymptotic expansion of a particular geodesic allows to recover the 
full picture describing the geodesic: that is the number of diffractions, the lengths 
of the legs between two diffractions, the diffraction angles and the angles of the 
cone at which the diffractions occur. We will tackle some of these questions in the 
second and third parts of this series. 

0.5. Organisation of the paper. In Section 1 we will recall the definition of sin¬ 
gular Fourier Integral Operators as defined in [MU79]. We will first study the case 
of two intersecting Lagrangians. We will give the oscillatory integral representation 
using a phase function that depends on an extra parameter s. We will then give 
the generalization to a system of four intersecting Lagrangians. 

In Section 2 we will study wave propagation on a cone of total angle Att. The 
first reason why we study this particular cone is that it is the simplest case in which 
we can implement our method of ‘moving the conical point’ that leads to our new 
expression for the wave propagator. The fact that the wave propagator belongs to 
the Melrose-Uhlmann class can then be directly read off from this expression. It 
is also worth remarking that, in this case the extra parameter s has a geometric 
meaning since it represents the amount of which the conical point has moved. 

The second reason why we can first study the cone of angle An is that the most 
singular part of the wave propagator near E actually does not depend on its angle. 
This can be seen using the construction of the wave kernel made by Friedlander in 
[FriSl]. We will recall this fact in Section 3 and then proceed to prove Theorem 

o~ 

In Section 4 we will use Theorem 0.3 to compute the wave propagator when 
microlocalized near some particular kind of geodesics. We will focus on the case of a 
geodesic with two geometric diffractions for which a desciption of the microlocalized 
propagator is not already available in the literature. 

In Section 5 we will end this paper by computing the leading contribution to the 
wave-trace of an isolated periodic orbit with two geometric diffractions. 

1. Intersecting Lagrangian distributions 

The class of distributions central to our study of the wave propagators on C'a is 
that of intersecting Lagrangian distributions, introduced by Melrose and Uhlmann 
[MU79]. These are distributions whose singularities (in terms of wavefront set) lie 
along a pair of conic Lagrangian submanifolds (Ao,Ai) of the cotangent bundle. 
Here, Ai is a manifold with boundary, and Aq and Ai intersect cleanly at 9Ai. 
In particular, the intersection is codimension 1 in both Lagrangians. These dis¬ 
tributions were introduced to construct fundamental solutions to operators of real 
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principal type. An analogous class of distributions associated to four intersecting 
Lagrangian submanifolds, also introduced in [MU79], will show up in our study of 
the wave kernel on a ESCS after two diffractions—see Section 1.3 . 

1.1. Model Lagrangian submanifolds. Let A be a manifold, and let (Aq, Ai) be 
a pair of conic Lagrangian submanifolds of T*X \ {0} with the geometry described 
above: Ai is a manifold with boundary, and Ag and Ai intersect cleanly at 9Ai. 
Moreover, let g S Ag n Ai be a point in the intersection. Melrose and Uhlmann 
showed that there is a normal form for this geometry. Indeed, let (Ag, Ai) be the 
model Lagrangian submanifolds in T*IR”: 

Ao = N*{0} = {{x,O--x = 0,} 

Ai = N*{x' = 0,a:i > 0} = {(a:,^) : x' = 0,^i = 0,xi > 0}. 

Here we decompose x = (xi,x'), where x' = (x 2 ,...,x„); similarly, ^ 

Choose any point q G Ag H Ai. Then Melrose and Uhlmann showed that there is 
a homogeneous sympectic map from a conic neighbourhood of g to a conic neigh¬ 
bourhood of q, such that A^ gets mapped to A^. To define intersecting Lagrangian 
distributions, they first defined them in the model situation. We recall this defini¬ 
tion. 

Definition 1.1 (Melrose-Uhlmann). An intersecting Lagrangian distribution of 
order m associated to the model pair (Ag, Ai) is a distributional half-density given 
by an oscillatory integral of the form 

( 1 . 2 ) (27r)-”-5 J 

where a is smooth, compactly supported in x and s, and a symbol of order m+ 5 — f 
in The space of such distributions is denoted /™(A; Ag, Ai). 

It is shown in [MU79] that elements of /"*(A;Ag,Ai) are Lagrangian distribu¬ 
tions of order m on Ai when microlocalized away from Ag, and Lagrangian dis¬ 
tributions of order m — ^ on Ag when microlocalized away from Ai. Also, they 
showed that the space /"*(A;Ag,Ai) is invariant under the action of Fourier in¬ 
tegral operators that fix Ag and Ai. Consequently, one can define intersecting 
Lagrangian distributions for a general pair (Ag,Ai) to be the image of the model 
space /“(A; Ag, Ai) under an FIO mapping A^ to A^. The precise definition is as 
follows: 

Definition 1.2. Let (Ag, Ai) be a pair of intersecting conic Lagrangian distribu¬ 
tions in T*X \ {0} with the geometry described above. The space /’"(A; Ag, Ai) 
consists of those distributional half-densities u that can be written as a locally finite 
sum 

U = Mg -I- Ml -I- ^ Fi{Vi) + Moo, 

i 

where Mg S /’"“^(XjAg), mi S /'"(A;Ai \ Ag), Vi G /“(A; Ag, Ai), Fi are FIOs 
mapping (Ag, Ai) to (Ag, Ai), and Uoo is C°°. 

In what follows, we will often omit the space ‘A’ from the notation for these 
distributions, i.e., we will write /™(Ag, Ai) in the place of /“(A; Ag, Ai). 
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1.2. Parametrization of intersecting Lagrangian submanifolds. Over the 
course of this paper, we will construct the fundamental solution of the wave kernel 
on a two-dimensional cone directly; we will want to be able to identify it as an 
intersecting Lagrangian distribution. To do this, we need a direct definition of 
intersecting Lagrangian distribution in terms of a phase function parametrizing a 
given pair (Aq, Ai) in place of the indirect Definition 1.2. 

Definition 1.3. Let (Ao,Ai) be a pair of intersecting Lagrangian submanifolds, 
and let q G ApnAi be a point in the intersection. A local parametrization of (Aq, Ai) 
near g is a function (j){x, 9, s), defined in neighbourhood of {xq, Oq, 0) C Ai xK'' xK^o 
such that 

• do, 0 ) = 0 , and q = {xq, d^(j>{xo, Oq, 0 )); 

• the differentials 


d(j) 


and dx,e 


d(j) 

ds 


(1.3) 


(1.4) 


in the {x,0) directions are linearly independent at (ccoj^ojO); 
the map 

Co {(x, 9) : dg(j){x, 9, 0) = 0} ^ {{x, d,^{x, 9, 0))} C T*X 

is a local diffeomorphism from Co onto a neighbourhood of g in Aq; 
the map 

Cl {{x, 9, s) : dg^s4>{x, 9,s) =0, s > 0} >-)■ {(a:, dx(t){x, 9, s))} C T*X 
is a local diffeomorphism from Ci onto a neighbourhood of g in Ai. 


Let us make some remarks about the definition above. The second condition 
ensures that the sets Cq is a smooth submanifold of Ai x of dimension n = dim Ai, 
and Cl is a smooth submanifold of Ai x x of dimension n transverse to 
{s = 0}. This makes it possible to speak of diffeomorphisms from Ci to A^ as in 
the third and fourth conditions. The first condition simply says that the base point 
(xo, 9o, 0 ) corresponds to the base point q. 

Proposition 1.4. (i) Let (Ao,Ai) C T*X be a pair of intersecting Lagrangian 
submanifolds, and let q be a point in the intersection. Then there exists a local 
parametrization o/(Ao,Ai) near q. 

(ii) Let (j), defined in a neighbourhood U of (xo,6>o,0) G Ai x x R>o, be a 
local parametrization of (Aq, Ai) near q. Let a(x, 9, s) be a classical symbol of order 
m — ^ + + j in the 9 variables which is compactly supported in U. Then the 

oscillatory integral 


(1.5) (2t)- 

is in /’"(Aq, Ai). 


/b'= Jo 


e^<l>C,o,^)a{x,9,s)ds d9\dx\i 


Proof, (i) By [MU79], there is a homogeneous canonical transformation y defined 
in a neighbourhood of 9 G Aq H Ai mapping Aq to Aq and Ai to Ai, and sending q 
to q. Let '^{x,y,9) be a phase function parametrizing the graph of this canonical 
transformation. Consider the sum of the phase functions 

'^{x,y,9) + yT]-r]iS, 
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where the second phase function is the standard parametrization of the model 
Lagrangian pair. Following [H6r71, p. 175], we define a new variable 

Y=\9\y. 


We then write this sum of the phase functions in terms of Y. That is, we define 


(j){x,Y,0,f],s) = ( a;, 1^,0 




Notice that 4> is homogeneous of degree 1 in the variables {Y,9,ri). We claim that 
(/) is a nondegenerate local parametrization of (Aq, Ai) near q. 

Let (?/Oj0 O)O) be the point corresponding to q and (xo,yo,0o) be the point cor¬ 
responding to (q, q) in the graph of x- Then dg^Y,r]4‘ = 0 and s = 0 implies that 
de'i'{xo,yo,do) = 0, yo = 0, d^^'(a:o, 0, 0o) = -0 and d2;4'(a;o, 0,0o) = a(0 , y) = g, 
so the first condition in Definition 1.3 is satisfied. 

We next check that the second condition is satisfied, i.e., that (f) is nondegenerate. 
To do this, we claim that the differentials 


and d^J^ 


09 J 


Oyi 


are linearly independent at {xo,yo,9o). This is a consequence of the fact that ik 
parametrizes A^, the (twisted) graph of the canonical transformation x^ which 
implies that the functions yi and dy.'^ are coordinates on A^. Using the diffeomor- 
phism between 

= {(a;,y,0) : dg^ = 0} 

and A^, we see that yi and are coordinates on C^s,. This implies that 


Vi, 


OVj’ 


and 


54' 


have linearly independent differentials at (xo,yoj^o)- Equivalently we can say that 


54- 

Oyj 


and d^g 


54- 


are linearly independent at (xq, j/Oi ^o)- This in turn is equivalent to the statement 
that 

dcj )' 


( 1 . 6 ) 


d(j) 


OY, 


and dxx 


are linearly independent at {xo,Yo,9o), 


where Yq = yol^ol- Now, from the explicit form of (j) it is evident that 


(1.7) 



and dY,Ti 



are linearly independent at (a;o,Voi^o)- 


Putting (1.6) and (1.7) together we find that ^ is a nondegenerate phase function, 
i.e., it satisfies the second point in Definition 1.3 . 

To check the third point, consider a point {x,Y,9,r],0) where dYfi^r]4> = 0 and 
s = 0. This implies that 


(1.8) dg'^{x,y,9) = 0, dy{y-r]) = 0, and dy'i'{x,y,9) + dy{y-r]) = 0. 
Using the fact that 4' parametrizes the twisted graph of X) this implies that 

(1.9) y = 0, dy^ = -ri, and (x, ^4') = x(y,-dy^')- 



WAVE PROPAGATION ON AN ESCS. I 


11 


Thus, dY,e,rj(t> = 0 implies that the Lagrangian parametrized is 

= {x(0,??)}. 

As {x,Y,9,r]) range over a neighbourhood of (tq, Tg, 6*o, lyo), the point (0,17) ranges 
over a neighbourhood oi q € Ag, and therefore x(0, rj) ranges over a neighbourhood 
oi q G Ag. This verifies the third condition in the Definition. Exactly the same 
reasoning shows that the fourth condition in the Definition is also satisfied. This 
completes the proof of part (i) of the Lemma. 

(ii) Choose an FIO F associated to the canonical relation x as above, and which 
is microlocally invertible at (q, q). Let F~^ denote a microlocal inverse to F. Write 
F~^ with respect to a phase function S{y,x,uj). Then the phase function 


$ = S{y,x,w) + (j){x,0,s) 

parametrizes the model pair (Ag, Ai) (after we homogenize the x variable by chang¬ 
ing to the variable X = a;|w|, as we did in the proof of part (i)). The proof is the 
same as in part (i), so we omit it. It then suffices to show that an oscillatory integral 
with phase function $, 

( 1 . 10 ) J F'^^y’^’‘^’^’^^a{y,X,uj,9,s)dsdXdeduj 

gives an element of /’"(Ag,Ai), since the original oscillatory integral is, modulo 
functions, the image of (1.10) by the Fourier integral operator F, which by 
definition maps /"“(Ag, Ai) to /™(Ag,Ai). Thus, we have reduced to the case that 
the intersecting pair (Ag, Ai) is the model pair (Ag, Ai). 

We now simplify our notation, and assume that ^{y,9,s) is a nondegenerate 
phase function parametrizing (Ag, Ai) locally near q, with (yg, 0g) corresponding to 
the point q. Here 9 G with k > n. We want to show that 

(1.11) « = for a S + t(x x K;^g;M'=) 

is in the space /™(Ag,Ai). Essentially this proof follows that of Proposition 3.2 

def 

in [MU79]. We first note that $g(y,0) = <I>(y, 0,0) parametrizes Ag. We have 
by [H6r71, (3.2.12)] that the rank of dgg^{yo,9o) is k — n. By rotating in the 9 
variables we can arrange that 9 — (9\ 9") with dim0' = n, dim0" = k — n and such 
that dg„g„^{yo,do) is nondegenerate. Integrating in the 9" variables and applying 
the stationary phase expansion, as in [H6r71, p. 142], we find that the result takes 
the form 


(1.12) ^ = /y for dGS 

where 9"{y^ 9', s) is the critical point, determined by the equation 

de"(t>{y,0' ,e" ,s) = 0; 




this varies smoothly with {y,9',s) near (yg,0g,O) thanks to the implicit function 
theorem and the nondegeneracy of dg„g„^{yo,9o) near (yg,0g,O). Then the phase 

function $g(y,0',O) $(y, 0', 0"(y, 0', 0), 0) parametrizes Ag. Moreover, it has 

the same number of hbre variables as the standard phase function y ■ y, and its 
fibre Hessian, dg'g/^Q has the same signature (namely, zero) as the fibre Hessian 
of y • 77. By Hormander’s equivalence of phase functions, [H6r71, (3.2.12)], there is 
a coordinate transformation 77 = r]{y, 9') mapping $g to y • 77 in a neighbourhood 
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of (2/O)^o>0)' Employing this change of variables, we are reduced to the case that 
$(y, 0', s) has the form y ■ 9' + 0(s). We can now follow the proof of Proposition 
3.2 in [MU79] from Equation (3.7) of [MU79] to the conclusion, which completes 
the proof of the Lemma. □ 


We next want to identify the symbols at Ag and Ai directly from the oscillatory 
integral expression (1.5). Recall that the symbol on each A^ is half-density taking 
values in the Maslov bundle. For our purposes, it is enough to do this when our 
Lagrangians Ag and Ai are conormal bundles. In this case, the Maslov bundle is 
canonically trivial, which means that we may regard the symbol as being simply a 
half-density. In the following theorem, we identify functions on A^ and Ci, where 
Ci is given by (1.3), (1.4). We let A = (Ai,..., A„) be local coordinates on Ci, or 
equivalently on the Lagrangian Ai. Similarly, we let A be local coordinates on Cg. 
Notice that we could choose A, A to be of the form A = (A',s) and A = [X'.dscj)) 
where A' are coordinates on Cg fl Ci. 


Proposition 1.5. Suppose now that Ag and Ai are both the conormal bundle of 
codimension one submanifolds Mq and Mi. Then 
(i) The symbol of (1.5) at Ai is given by 


(1.13) 


a{x,9,s)\ci 


d{X,(j)e,(ps) 

d{x, 0, s) 


\dX\i 


where a is the signature of the Hessian (j)”g in the {0, s) variables, 

(ii) The symbol of (1.5) at Ag is given by 


(1.14) 


(27r)- 


ia{x, 0, 0) 

4>s(.x,0,Q) 


1 

d{X,(j)e) 

ICo 

d{x,0) 


1 

2 


|dA|5 


where a is the signature of the Hessian (jfgg at s = 0. 


Remark 1.6. We remark that cr and a' are constant, as follows from [H6r71, 
(3.2.10)] by comparing with the standard parametrization of a conormal bundle 
with linear phase function. 


This proposition follows directly from [H6r71, Section 3]. 

1.3. Four intersecting Lagrangians. The wave kernel after two diffractions is 
associated to four different Lagrangian submanifolds: the direct front, one front 
from a diffraction with each cone point, and a fourth front from diffractions with 
both cone points. We shall show that the wave kernel in this case is contained in 
the Melrose-Uhlmann calculus of distributions associated to four Lagrangian distri¬ 
butions described in [MU79, Sections 7-10]. We now recall some of this material, 
starting with the definition of a system of intersecting Lagrangian submanifolds. 


Definition 1.7. A system of four intersecting conic Lagrangian submanifolds of 
T*X is a quadruple A = (Ag, Ai, A2, Ag) of Lagrangian submanifolds, where Ai 
and A2 are manifolds with boundary and A2 is a manifold with codimension two 
corner, with the following properties: 

• (Ag,Ai) and (Ag,A2) are intersecting pairs in the sense of the previous 
subsection; 

• Ai n A2 = i9Ai n 9A2 = Ag n Ag = cAg, where cAg denotes the codimension 
2 corner of Ag; 

• The two boundary hypersurfaces of Ag are Ag n Ag and Ag n A2. 
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For example, the following is a system of intersecting Lagrangian submanifolds: 

Definition 1.8. Suppose n ^ 3. For j = 0, define A = (Ag, Ai, A2, A3) to 

be the following quadruple of Lagrangian submanifolds of T*IR": 

o,a = 0} 

Xn = 0,^2 = 0 } 

= 0,^1= 6 = 0}. 

Locally, an intersecting system as in Definition 1.7 may be realized as follows. 
Let Aq be a Lagrangian submanifold, and let pi, p 2 be two functions on T*X such 
the Hamilton vector fields L/pi, Hp^ are linearly independent, transverse to Ag, 
and commute with each other. Then we define A^, f = 1, 2 to be the flowout from 
Ag n {pi = 0} by Hp-, and A3 to be the flowout from Ag n {pi = P 2 = 0} by 
the flowout of both Hamilton vector fields. For example, the model system is of 
this form, where pi = 6 P2 = 6- It turns out that, locally, all intersecting 
systems arise in this way. As a consequence, every system of four intersecting 
Lagrangian submanifolds is the image of a model system under a homogeneous 
canonical transformation. We now define the model system. That is, one could 
alternatively define an intersecting system by the requirement that, locally, it is the 
same of the model system under a homogeneous canonical transformation. 

We next define the space of Lagrangian distributions associated to the model 
intersecting system A given by (1.15). 

Definition 1.9 ([MU79, Definition 8.1]). The space A) consists of those 

distributional half-densities that can be expressed in the form 

n pOO pOO 

(1.16) (2 t)-”-M / / e^^^-^-^^^^-^^^^^a{x,^,si,S2)dsids2d^\dx\^ 

J Jo Jo 

where a is smooth and compactly supported in (x, si, S2) and is a symbol of order 
m -I- 1 — ^ in the ^variables. 

It is not hard to check that if rt G /™(IR"; A) then the wavefront set is of u is 
contained in IJj^g A^, and if g G A^ is not contained in Aj for j 7^ i, then m is a 
Lagrangian distribution associated to A^ microlocally near q, of order to if f = 2, 
TO—Iiff = lor2 and to— 1 if i = 0. We can also observe that if u is microsupported 
near A^ H A^ , i < j, and away from the other A^, then it is an intersecting pair 
of order to — | associated to (Ai,Aj) for {i,j) = (0,1) or (0,2), or of order to for 
(i, j) = (1,3) or (2,3). 

It is shown in [MU79] that the model space /^(K"; A) is invariant under FlOs 
that map each A^ to itself. As a consequence, we can define intersecting Lagrangian 
distributions associated to a general intersecting system A = (Ag, Ai, A2, A3). 

Definition 1.10 ([MU79, Definition 8.7]). Let A be an intersecting system of 
homogeneous Lagrangian submanifolds of T*X. The space /’"(A; A) consists of 
those distributional half-densities u that can be written as a locally finite sum 

u = uoi+ U 02 -I- mi 3 + U23 + X! 


Ao = {(a:,6 : x = 0} 

Ai = {(a :,6 : xi > 0,2:2 = •••= a;„ = 
A2 = {(a:,6 : 2:2 > 0,a:i = 2:3 = • • • = 
A3 = {(a:, 6 : xi > 0,2:2 > 0,2:3 = ... 
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where Uij S /"* for (i,j) = (0,1) or (0,2), Uij £ Ai, Aj) for 

ihj) = (Ij 3) or (2, 3), Fi are FIOs mapping the model intersecting system A to A, 
and Vi e A). 


As before, we will often omit the space ‘A’ from the notation for these spaces of 
distributions. 

We will find it useful to have a definition of /'"(A) defined directly in terms of 
phase functions. To this end we give an analogue of Proposition 1.4 in the setting 
of intersecting systems. We first need a definition of a phase function parametrizing 
an intersecting system A, locally near a point q G A. Notice that either q is in only 
one of the A^; or in one of the four-fold intersections Aq H Ai, Ag H A 2 , Ai n A 3 , 
or A 2 n A 3 , and disjoint from the other two; or in the 4-fold intersection 
Since these four pairs (A^, A_,) form intersecting pairs of Lagrangian submanifolds 
in the sense of the previous subsection, the only case in which we have not already 
defined a local parametrization is in the case that q £ nloA.. 


Definition 1.11. Let A = (Ag, Ai, A 2 , A 3 ) be a system of intersecting Lagrangian 
submanifolds, and choose a point q £ n^=c Ai in their intersection. We say that 4> 
is a local parametrization of A near q if it is a function (j>{x, 9, si, S 2 )) defined in a 
neighbourhood of (xg, 0o, 0, 0) C M x (K^ \ {0}) X ^ homogeneous of 

degree 1 in ^ such that 

• de,si,s 2 'P{xo, Oo, 0,0) = 0, and q = {xq, d^cj^ixo, 6»g, 0,0)); 

• the differentials 


(1.17) 



and dr 


d<j) 

ds2 


in the {x, 0 ) directions are linearly independent at (a;g, 0g, 0 , 0 ); 
• the map 


(1.18) 


Co ='{(ai,0) 


de<l)ix, 9, 0,0) = 0} hA {(a;, d^<))(a;, 9, 0, 0))} C T*X 


is a local diffeomorphism from Cq onto a neighbourhood of g in Ag; 

• the map 

Cl = {(a;, 6 », si) : de,si4>ix,9,si,0) = 0, si > 0} >-)■ {{x,dr(l){x,9,si,0))} C T*X 

is a local diffeomorphism from Ci onto a neighbourhood of g in Ai; 

• the map 

C 2 {(a;, 9, S 2 ) : de ,82 4>{x^ 0, S 2 ) = 0, S 2 > 0} >-)■ {(a:, dr4>(x, 9, 0, S 2 ))} C T*X 

is a local diffeomorphism from C 2 onto a neighbourhood of <7 in A 2 ; 

• the map 


C 3 = {(a;, 6 », si, S 2 ) : dg^si,s24‘{x, 9, 81 , 82 ) =0, si > 0, S 2 > 0} 

1-4 {{x,dr(l)ix,9,Si,82))} 

is a local diffeomorphism from C 3 onto a neighbourhood of <7 in A 3 . 

Proposition 1.12. (i) Let A = (Ag, Ai, A 2 , A 3 ) C T*X be a 8ystem of intersecting 
Lagrangian submanifolds, and let q be a point in the intersection. Then there exists 
a local parametrization of A near q. 
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(ii) Let (j), defined in a neighbourhood U of (xq, Oj 0) be a local parametrization 
of A near q. Let a{x, 6, si, S 2 ) be a classical symbol of order m—| + 1 + ^ in the 
0-variables which is compactly supported in U. Then the oscillatory integral 



is in 7™(A). 

Proof. The Proposition is proved in the same way as Proposition 1.4. □ 

Remark 1.13. For a given phase function 9, si, S 2 ) to parametrize some sys¬ 
tem of four intersecting Lagrangian submanifolds, locally near (ccq, Oj 0), it is 
sufficient that it satisfies dg^si,s24>{xo:(^o,0,0) = 0 and condition (1.17). Then the 
sets Ai, i = 0 .. .3, defined as the image of Ci in Definition 1.11, are automatically 
Lagrangian submanifolds satisfying the geometric conditions to form a system in 
the sense of Definition 1.7. 


We next write down an expression for the symbol of the oscillatory integral (1.19) 
at Aq. As in the previous section, we restrict to conormal bundles, in which case 
the Maslov bundle is canonically trivial. We write A = (Ai,..., A„) for coordinates 
on Co which we identify with Aq via (1.18). 


Proposition 1.14. Using the notation of (1.19), suppose now that Aq is the conor¬ 
mal bundle of a codimension one submanifold Mq and Mi. Then the symbol of 
(1.19) at Aq is given by 


( 1 . 20 ) 


- ( 271 -)’ 


a{x, 9, 0 , 0 ) 




Co 


d{\fe) 
d{x, 9) 


_(/)si(a;,6>, 0, 0)(^s2(a;, 6», 0, 0)_ 

where a is the signature of the Hessian fgg at Si = S 2 = 0, and 
local coordinates on Cq. 


2. The microlocal structure of the wave propagator on C4,r 


We now specialize to the cone C 4 ,r, where we will carry out the actual analysis 
of the sine propagator near the singular set. Let us first pause for a moment to 
highlight some features of the cone First, and perhaps most important, the 
interior C|,r is equivalent to the double cover of the punctured plane R^\{(0, 0)}. As 


def 

a result, the Schwartz kernel E = K. 


i(tVA) 


vW 


in this setting (cf. [CT82b, p. 448-9]): 

( 2 . 1 a) E{t,ri,9i-r2, 


has a particularly simple description 


2 ) = 0 


when 0 < t < dist(ri, 6 <i; r 2 , 6 * 2 ); 

(2.1b) E{t, ri, 6 »i; r 2 , 6 * 2 ) = ^ “ (^1 + '’2 “ 2 rir 2 cos( 6 »i - 6 » 2 ))] ”" 

when dist(ri, di; r 2 , ^ 2 ) < ^ < -I- r 2 ; and 

( 2 . 1 c) F;(t,ri, 6 ii;r 2 , 6 » 2 ) = ^ - (r-i-f r 2 - 2 rir 2 COs( 6 ii - 02 ))] ", 

when t > ri -\- r 2 . In particular, the jump discontinuity across the diffractive front 
{t = Ti -I- r 2 } is readily apparent on C 4 .n..^ Second, a seemingly incidental fact 


^Note that [CT82b, eq. (4.7)] contains a sign error that we have corrected here. 
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that will be important as we continue is that constant vector fields are well-defined 
on ( 74 ^ (and indeed any cone with cone angle an integral multiple of 27r, i.e., the 
finite-sheeted covering spaces of the punctured plane). 

2.1. The ‘moving conical point’ method. Our technique for determining the 
structure of the wave kernel is the ‘moving conical point’ method. Given two 
points qj and in and a positive time t* , we want to determine E{t, qi, 92 ) 
for {t, qi, q 2 ) in a neighbourhood of (t*, q*, qj)- To do this, we imagine that we can 
move the conical point (that is, the place where the two copies of are ramified) 
along a straight line, in a direction such that moves it ‘in between’ q* and q^, 
and then far away (i.e., at a distance S much larger than t*). This means that 
the angle between q* and q^ tends to 27r, so the distance between them will be 
2S + 0(1) S> t* . Then, by finite propagation speed, after the cone point is so 
shifted, the wave kernel at (f*, q*, q|) will vanish. We then express the kernel using 
the fundamental theorem of calculus: 

/■■S' 

E{t,qi,q2) = - J —E‘^{t,qi,q2)ds, 

where E^{t,qi,q 2 ) is the wave kernel where the cone point has been shifted a 
distance s in our chosen direction. Thus, if we can understand the derivative of E'^ 
with respect to s, then we can compute E = E'^. The reason we can expect the 
derivative -^E^ to be simpler than E^ itself is that the singularity at the direct front 
is independent of s, so ^£ 1 ® should be associated purely to diffractive behaviour. 
The rest of this section is devoted to implementing this method. 

To do this in a rigorous manner, rather than moving the cone point, we instead 
translate the points qi and q 2 on the cone (in the opposite direction — see Fig¬ 
ure 2.1) using the flow of a constant vector field X G which we choose in 

a direction such that the two half-lines parallel to X through qi and q 2 pass on 
different sides of the cone point; in particular, neither meets the cone point. 

We set (p^ = (fix to be the associated flow, the group of tocat^ diffeomorphisms 
given by time-s translation along X. Using we assemble the kernel spacetime flow 
for X, which is the group of locally-defined diffeomorphisms on Rf x x 
given by 

^‘‘{t,qi-,q 2 ) {t,(p"{qi);‘p%q 2 )) = {t,qi + sX;q 2 + sX). 

Consider the distribution 

( 2 . 2 ) Es=\ds[{^^yE], 

where y = x(qi; q 2 ) G C°°(C'|^ x Cl^) is a smooth function that vanishes near the 
cone point. Its role is to ensure that (<i)®)*£l is well defined on the support of y; 
that is, X must be chosen so that it vanishes in the set obtained by translating a 
small ball centered at p in the X direction, and is identically 1 in a neighbourhood 
of the set {(qi -|- sX, q 2 -I- sX) : s G K, (xi,a; 2 ) G U}, where 17 is a suitably small 
neighbourhood of {q*,qy. Then for (qi,q 2 ) G U we have 

(2.3) s, = X • a, [{^y*E] = ds [($^)*£;]. 


^The time interval for which cp^ is defined depends on the starting point; in particular, the points 
along the reverse fiowout of p can only be evolved forward for finite time—until they reach p. 
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Figure 2.1. Moving the conical point. Shown are the singular 
support of E^{t,qi;q 2 ) = E{t,qi{s)]q 2 {s)) (solid circles) and the 
singular support of E{t, qi; 92 ) (dotted circles) as the moving con¬ 
ical point p{s) travels along the flow of —X; the branch cut is 
depicted as the red dashed line. 


Set 

(2.4) Ts{t,qi;q2 ) x{QuQ2 ) ■ (i, Qi; 92 ); 

this is the precise version of the quantity -^E^ in the heuristic discussion above. 
Thus, we have 

(2.5) E(t,qi;q2) = - f T ,{t, qi, q2) ds, ( 91 , 92 ) et/, 

Jo 

provided that {^^)*E{t, 9i; 92 ) = 0 as discussed above. 

When s = 0, we calculate that 

(2.6) To(t,9i;92) = XiF;(t, 91192 ) -fX2Fl(t,9i; 92 ) 

with Xj denoting X acting in the 9 j-variable, and for general s we have 

(2.7) 9i; 92 ) = (4'®)*To(<, 9i; 92 ), ( 91 , 92 ) £ U, 
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since the vector field X is constant. Pairing Tq with a test function ip G 
in the q 2 "Variable, we then integrate by parts to obtain 

{To,iP)^^ = {X^E,iP)^^ + (X2E,iP)^^ 

= {X^E,iP)^^-{E,XiP)^^ 

= (X o W{t))iP - (W(t) o X)iP 
= [X,W{t)]ip. 

Thus, To is the Schwartz kernel of the commutator [X, W(t)] of the constant vector 
field with the sine propagator. Note this distribution is everywhere well-defined. 

A quick computation now yields the operator identity 

□ o [X,W(t)] = [X, A] o W(t), 

and hence Duhamel’s principle implies^ 

(2.8) [X, W{t)] = - f W{t -s)o [X, A] o W(s) ds, 

Js=0 

where we recall the Schwartz kernel of these operators is Tq. Using (2.8), we will 
show that Tq is a multiple oi S{t — ri — r 2 ), hence a purely diffractive Lagrangian 
distribution. First, we must understand better the commutator [X, A]. This is the 
aim of the next subsection. 


2.2. Distributions supported at the cone point and commutators. To make 
full use of the expression (2.8), we need to understand explicitly the Schwartz kernel 
of the commutator [X, A]. This requires a brief detour through the spectral theory 
of the Laplacian on C 4 ,r, and, in particular, a discussion of the failure of essential 
self-adjointness of the Laplace-Beltrami operator Ag on 

Let iL®(C' 47 r) denote the usual Sobolev spaces on € 4 -,^, defined as 

(2.9) =' {u e L\C4.) : Diff'=(C' 4 .) • u G ^^(Q.)} 

for integers k G and extended to all real orders by duality and interpolation. 
An exercise (essentially the same as a more standard calculation on where 
the same result holds; cf. Chapter 1.5 of [AGHKH05]) shows that the closure of 
in the graph norm for Ag, 

|k||A/= iklU^+llAgUlU., 


i.e., the domain of the closure Ag of Ag, is 


(2.10) D =^Dom(Ag) = {u G : u{p) = O} . 

Thus, if p G Cp°{C 4 Tr) is any bump function satisfying p = 1 for r ^ 1 and p = 0 
for r ^ 2, then this shows 

H‘^{C4t,) = 2) © Spanc{p}. 


We show in Lemma A.l that the domain of the adjoint of this operator is 
25* Dom(Ag*) = 2) © Spanpl p, plog(r), pr^ exp 




,pr 2 exp 




^There is a minus sign in the formula because our operator □ is — A = —d^ — A, while the 
usual Duhamel formula is written for an operator with a positive sign in front of . 
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The choice of a self-adjoint extension of Ag is then the suitable choice of a half¬ 
dimensional subspace of D*/D (cf. [RS75] for more details on self-adjoint exten¬ 
sions) . 

In our analysis, we have elected to work with the Friedrichs extension A 
of the Laplacian, the unique self-adjoint extension whose domain contains the form 
domain (which in our setting is H^{CiT^)). We define the spaces Vg to be the 
domains of real powers of this operator: 

(2.11) =^Dom(A5) . 

For s > 1, these spaces are strictly larger than the Sobolev spaces In 

particular, I ?2 is the Friedrichs domain itself. 

To distinguish the elements of I ?2 from those of 2), we must examine their be¬ 
havior at p. We do so in the following lemma, which we prove in Appendix A. 


Lemma 2.1. Fix a compactly supported, smooth, and radial cutoff p € C(?°(C' 4 ^) 
which is identically 1 near p. For any function u S I? 2 ) there exist constants a_i, 
oq, and oi in C and a distribution u S 2) such that 


( 2 . 12 ) 


= ( oo -I- a-i-sA" exp 


- 2 ^ 


-I- oi y/r exp 


2 ^ 


p{r) 


In particular, the function u — ao — O-i^/r exp[—| 0] — oi-^/r exp[| 0] vanishes at 


P- 


Remark 2.2. We see from Lemma 2.1 the system of strict inclusions 

® C C X>2 £ So¬ 


using this lemma, we see that the Friedrichs extension exactly corresponds to 
the choice of the functions 


Mr,0) = h 



and 


F+iir,0) 


def /— 



as the models for the admissible asymptotics at p. Given a function u in I? 2 , we 
define the distributions Lj for j = —1, 0, or +1 as 


(2.13) Lg{u) = ag 

in terms of the expansion (2.12). Note that the expansion (2.12) is independent of 
the choice of the cutoff p, for the difference of any two such cutoffs is compactly 
supported in and is thus in J). Hence, the distributions Lj are well-defined 
elements of J)_ 2 . Equivalently, we may define the Lfs using the angular spectral 
projectors 


(2.14) 


[n^u](r) 



u{r, 9) exp 



d0, 


and a straightforward computation shows that 
(2.15) ' ' ' ^ 


Lo{u) = . — lim [HomI (r) and L±i{u) = . — lim 

V 471 NO -'—''' 


[n±iu] (r) 


r-|o sfr 


Directly from the definition or from the above, we observe that L±i(u) = Lz^i(u). 

Corollary 2.3. Suppose L is a distribution in I ?_2 which is supported only at the 
cone point p. Then L is a linear combination o/L_i, Lq, and Li. 
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Proof. Suppose u is an element of T> 2 - By (2.12) we have 

L{u) = aoL{p{r)) + a_i L(<p_i(r, 6») • p{r)) + ai L{(pi{r,0) ■ p{r)) 
since v being an element of 2) implies L{v) vanishes. Therefore, 

L = L{p{r)) ■ Lo + 0) ■ p{r)) ■ L_i + L{pi{r, 0) ■ p{r)) ■ Li, 

showing that rt is a linear combination of Lq, T-i, and Li as claimed. □ 

Returning to the commutator [X, A], let us observe 

X:V2^Vi=H\C4^) 

since I ?2 is not contained in On the other hand, since is con¬ 
tained in I? 2 , we certainly have X : I ?2 —5" and hence, by duality, also 

X : L^(C' 4 ,r) —>■ 25-2- Therefore, for any u S V 2 , the commutator [X, A] is in I?_ 2 . 
On the other hand, if u is compactly supported in C 4 ,,., then the action of A on u 
is the same as the Euclidean Laplacian acting on u. Since the Euclidean Laplacian 
commutes with constant vector fields, this implies [X, A] m = 0. Therefore, the 
distributional support of [X, A] u for any u € I ?2 is at most the cone point p, and 
thus it fits into the framework of Corollary 2.3. 

Proposition 2.4. Let X = + X^d^ be a constant vector field on 

written in terms of the complex coordinate w = x + iy = re*®. Then for any 
distribution u G Vk for k ^ 2, we have 

(2.16) [X,A]m = -27r(A„Li(u) -Li -t AwT_i(m) -L.i) 

Proof. Consider the bilinear pairing 

([X, A]M,r), u,vGT>2- 

The discussion above shows that this is well defined for all u,v G T> 2 - It is clear 
that this pairing vanishes if either u 01 v lie in 2). So to compute the pairing, it 
suffices to consider u and v to be linear combinations of the functions Co = p(l'f^P)j 
C-i(w) p(|wp) and Ci{w) p(|wp). In fact, the pairing also vanishes if 

either u or v are Co since this is equal to a constant in a neighbourhood of the cone 
point, hence vanishes near the cone point after the application of either X or A. 
So we need only consider u and v equal to a combination of C±i- 

First, let X = For this X, consider the action of the operator [X, A] on a 
Fourier mode for a half-integer j. Since the Fourier modes are eigenfunctions 
of Agi , and since dw maps e*-^® to a multiple of the same property is true 

of [X, A]. It follows that the only nonzero combination with X = 9m is 

([9m, A]Cl,Cl)- 

Similarly, when X = 9m, the only nonzero combination occurs when u = u = C-i- 
In view of these considerations, to establish (2.16), it suffices to show that 

([5m,A]Ci,Ci) =-27r, 

([^, ^]C-i, C-i) = —27r. 

In fact, as the cafcufations are similar, we only prove the first. 

Since we are using the bilinear pairing we have 

([9m, A]Cl,Cl) = -2(ACl,9mCl) = -2 / ACl9mCld^ 
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where dS denotes the Euclidean area element. Using Stokes formula we have: 
[ {dwCifdw = ij ACi5u,Ci 

J r—e J r'>s 

For e small enough we thus obtain 

[ A(:M,ds = 7T. 


/ r>£ 


The claim follows. 


□ 


2.3. The differentiated wave propagator on (74^. We now apply the formula 
(2.16) for the Schwartz kernel of the commutator [X, A] to the Duhamel formula 
(2.8) to compute the distribution Tq. Writing X in complex coordinates, i.e., 
X = Au, dw + dw, this yields 


(2.18) To(t, gi; 92 ) = 2^ / {x^ [W(t - s)Li] (gi) • [Li o W(s)] {q^) 

Js=0 

+ X^[W{t - s)L_i] (gi) • [L_1 o W(s)] (g2)} ds. 

In particular, this shows Tq is an integral superposition of tensor products of the 
distributions 

(2.19) 

obtained from evolving the distributions Lj under the sine flow W(t). (Note that 
the self-adjointness of W(t), and the fact that its kernel is real, implies W(t) Lj = 
Lj o W(t), so we only need to work with the evolved distributions £j{t).) Since 
the Lj's are supported only at the cone point p, we should expect the propagated 
distributions £j{t) to be spherical waves emanating out from p, i.e., they should be 
diffractive-type waves. As the next lemma shows, this is indeed the case. 


Lemma 2.5. Let t > 0. The distributions ii{t) and £-i{t) on are given 
explicitly by 


( 2 . 20 ) 


£±i{t) = 


iTT-y/r 


6(t — r) exp 




Proof of Lemma 2.5. It suffices to prove the lemma for £i (t) ; the statement for 
£-i{t) is similar and follows by complex conjugation. 

Recall the spectral projector form of the dehnition of Li (see (2.15) and (2.14)) 
and 


Li{u) = lim-—^ 
rio AiTyr 


u{r, 6) exp 



dO. 


To compute the action of W(f) on Li, we use Cheeger’s functional calculus on 
metric cones [CT82a]; this expresses E as the sum 


-E(t,ri,6»i;r2,6>2) 


1 

4?! 


^exp 

jei- 


2 


{0i 



sin(At) 


Jiii(Ar2) AdA 


A=0 


A 
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over the angular modes of A. Since Li vanishes except at the j = 1 mode in this 
sum, we have the following simple formula for the action of 


[ei{t)] (u) 


c 4 : 7 T pOO 




= lim 


rito (47r)2y^ 


sin(At) 


Ji (Ari) Ji (Ar 2 ) A dX 


0 LJA=0 


X exp 


— 6 »i 
2 


exp 


2(^1 -^ 2 ) 


u{r 2 ,02) dOi r 2 dr2d62 


Performing the 0i-integral, this simplifies to 


(m) = liin 


poo pAir 


sin(Af) 


ritO AlTydri Jr2=0 7^2=0 I A=o 

X exp 


Ji (Ari) Ji (Ar 2 ) A dX 




u(r2,02) F2 dr2d02- 


We now substitute the explicit formula Ji{z) = [^] ^ sin(z) into the above, giving 


£i{t) = lim 


nto 27r2^ 


sin(Af) 


sin(Ari) 

Ari 


sin(Ar 2 ) exp 




dX. 


By pairing with a test function and using dominated convergence, we see that this 
is equivalent (in the sense of distributions) to the expression 


hit) = 


27r2^ 

To conclude the proof, we observe 


sin(At) sin(Ar 2 ) exp 


— 02 
2 


dX. 


e-*A(t+r) ^ 




/A =0 


f X——00 


This implies, dropping the subscripts from the base variables and replacing the sine 
functions by their complex exponential definitions, that 


lW = - 


1 


( 2 . 21 ) 


J\=- 

1 


I giA(t+r) _ ^iX{t-r) 


'1 

i „ 

I exp 

--0 

2 


dX 




S{t — r) exp 


— 0 
2 


, for t > 0 . 


□ 


Remark 2.6. It is remarkable that, on the cone of angle 47 r, there are solutions 
to the wave equation, namely r~^6it — r)e^*®/^ obeying the sharp Huygen’s prin¬ 
ciple, that is, supported on the light cone itself. This can be confirmed by direct 
calculation, applying the wave operator to these distributions. 

We also remark that one can prove Lemma 2.5 without appealing to the Cheeger 
functional calculus: after verifying that the 7±(t) satisfy the wave equation, it only 
remains to check that limi_>o ^±i(t) = 0 and \iuit^Q(d/dt)£±i{t) = L±i. 

We conclude this subsection with the proof that Tq is a Lagrangian distribution 
associated to the diffractive Lagrangian relation 
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Proposition 2.7. Let t > 0, and suppose, as in the diseussion in Section 2.1, that 
X points in the direction 9. Then the distribution Tg = /C [ [X, W(t)] ] is given 
explicitly in polar coordinates by 

(2.22) To(t,qi;q2) = - ra) cos^ - 9 ^ . 

Proof. We begin by rewriting the equation (2.18) using the distributions £j{t): 


T'o(t,gi;q2) = 271 [ |x„ [4(t-s)](qi) • [£i(s)](q2) 

Js^O ^ 

+ X ^[£^ i{t - s )]{ qi ) ■ [^_i(s)] (q2)| fis. 
We break up the integral across the sum and consider the first summand: 


Tg (t,qi;q 2 ) [ [£iit-s)]{qi) ■ [£i{s)]{q 2 )ds. 

Js=0 

Substituting our expression (2.21) in for £i{t) and its conjugate for £-i{t), the above 
becomes 


(2.23) 


1 /■* 

——= / S{{t - s) - ri)S{s - r 2 )exp 
\/^ 1^2 Js=0 


leir'^y/nr^ y^^g 

1 


WTT^y/rfrf 


S{t — ri — ra) exp 


-^( 01 + 02 ) 
{9i+92) . 


ds 


Similarly, we have 
(2.24) T(f(t,qi;q 2 ) 


IbTT^^ynra 


S{t — ri — ra) exp 


2 (01+02) 


For the vector field with direction 9, we have Xyj = exp{i9) = X^. Adding (2.23) 
times Xyj to (2.24) times Xi^jj, and multiplying by 2?!, we obtain (2.22). □ 


2.4. The full wave propagator on C' 4 ,r- Having computed Tg(t, qi; < 72 ), we re¬ 
turn to (2.5) and compute the sine wave kernel E[t,qi,q 2 ) on C' 47 r. Since our 
primary interest is in the behaviour near a geometric diffractive geodesic, let us 
assume for a while that 9i is close to 0 and 02 is close to tt (so that the diffraction 
angle 0i — 02 is close to —tt). We then choose to move the conical point in the 
direction This amounts to putting 0 = — f in the previous formulas. 

Let rj(s),0j(s) be the distance and angle from the point qj to the cone point 
shifted by a distance s in the 0 = f direction, or equivalently, from the point qj{s), 
obtained from qj by shifting a distance s in the 0 = — § direction, to the (fixed) 
cone point. Notice that, in the limit s —^ 00 , the angle between qi(s) and 52 ( 5 ) 
approaches 2?!. In particular, the points will be distance ri(s) -I- r 2 (s) = 2s + 0(1) 
apart, in this limit. Thus, the condition that {^‘^)*E{t, qp, ^ 2 ) = 0 is valid for large 
s. Hence, we can write, using (2.23) and (2.5), 

(2.25) 

E{t,qi,q 2 ) = ^ J (ri(s)r2(s))"5(5(t-ri(s)-r 2 (s))sin^^i^^^^^^^ ds. 

This can be written 

(2.26) E{t,q,;q2)= [ 

Js>0 




a(t, qi, q 2 , s, oj) ds doj 
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with the following phase function and amplitude: 


(2.27) 


(j){qi,q 2 ,s,uj) = [\Jx\ + (yi - s)^ + \Jxl + (2/2 - s)^ - 

0i(s) + 6*2(3) 


a{t,qi,q2,s,uj) = ^ • (ri(s)r2(s)) ^ • sin 


Since the phase function is a nondegenerate phase function in the sense of Dehni- 
tion 1.3 , we hnd that the propagator is in the Melrose-Uhlmann class. 

This construction can actually be carried out as long as S (—f j f) and 62 G 
(|, ^). When 6*1 is in the same interval but 62 now belongs to (—^,—f) (thus 
containing the diffraction angle of +7r), the conical point have to be moved in the 
opposite direction 6 = This leads to a similar expression. Observe however 
that in that case the phase is now 

<P{qi,q2,s,uj) = [y/xl + (yi +s )2 + \Jxl + (j/2 + sY - 

In the remaining cases for which 6*2 belongs respectively to (—|, |) and (^, ^) 
the conical point can be moved in the 9 = n direction. It should be noted however 
that in this case the limit s —>■ 00 of £^®(t, Qi, 92 ) is not 0 but the free solution. 

In any case, it follows that E is an intersecting Lagrangian distribution in a 
neighbourhood of (6, <71; 52 )- Close to the diffraction angle —tt, we use the form of 
the phase (2.27) to determine the two Lagrangian submanifolds. First, when s = 0, 
it is clear that (j)\s=o parametrizes the Lagrangian N*{t = ri + r 2 } = Second, 
when s = 0, ()) is stationary with respect to s when the cone point lies on the straight 
line between (xi,yi + s) and {x 2 ,y 2 + s). In this case, the second derivative 
is nonzero, and we can eliminate the variable s by replacing it with its stationary 
value. In this case, the sum of distances xf + {yi + s)^ + y/a;| + (?/2 + s)^ is 
equal to the distance between {xi,yi + s) and ( 0 : 2 , 2/2 + s), which is the same as the 
distance between qi and 92 - So an equivalent phase function is (|qi — < 72 ] — 
and this parametrizes the conormal bundle of the direct front, A*^. 

This essentially proves 


Proposition 2.8. For each fixed 6 > 0, the sine propagator kernel E on is an 
intersecting Lagrangian distribution on x C1^ of order — !.• 

E{t) G rP x A°, A^) ; 

in particular, it has Lagrangian order —1 on A*^ \ A*^ and order — | on A*^ \ A^. 

2.5. The Cheeger-Taylor formula. It is instructive to compute the integral 
(2.25) explicitly, and confirm that we obtain the Cheeger-Taylor formulae for the 
wave kernel from Section 2.1. Let us consider the case in which 61 G (—f, f) and 

^2G(f,f) ~ 

Since the functions ri{s) take the form a/fq + (s — sq)^, they are convex func¬ 
tions of s. Therefore, as s ranges from 0 to 00 , the delta function 5(t — ri{s) — r 2 {s)) 
can be nonzero for at most two values of s. More precisely, if 6 < ri -f r 2 and the 
angle between qi and q 2 is greater than tt, then there are no values of s for which 
t = ri(s) -l-r 2 (s), since in this case, both ri(s) and r 2 {s) are increasing in s. On the 
other hand, suppose that t < ri-\-r2 and the angle between Xi and X2 is less than tt. 
We might as well assume that t > d{xi, X 2 ), since otherwise the wave kernel is zero 
due to finite speed of propagation. In this case, ri(s)-|-r 2 (s) decreases until the cone 
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point lies directly between xi and X 2 j when we have ri(s) + r 2 (s) = d{xi,X 2 ) < t, 
and then increases to infinity. It follows that in this case there are two values of s 
for which t = ri(s) + r 2 {s). The final case is t > ri + r 2 . In this case, regardless of 
whether ri(s) + r 2 (s) initially increases or decreases, there is always one value of s 
for which t = ri(s) + r 2 {s). 

For each value of s satisfying t = ri(s) + T 2 (s), we calculate the contribution to 
the integral (2.25). This is given by 

(2-28) ^(ri(s)r2(s))"5|r;(s) + r2(s)|”Sin^Mfi±Mfi 


= ^(ri(s)r 2 (s)) 2 |sin( 6 »i(s)) + sin( 6 » 2 (s))| ^ sin 

47r 


= ^ [(?'l(s)’'2(s)]’ 


0l(s) + 6>2 (s) 


sin 




(s) + ^2(s) 

2 


sin( 0 i(s)) + sin( 02 (s)) 


since by choice ^d«)+^ 2 (a) ^ ( 0 , 7 r). Using the addition formula for sin0i + sin 02 we 


obtain that the contribution can be written 


—(ri(s)r2(s)) 5 


COS 


-1 


and we want to prove that this coincides with 

^ (f - {rl + r 2 - 2 rir 2 cos( 0 i - O 2 )) 


whenever the moved conical point p{s) lies in between qi and < 72 . This implies that 
t = ri{s) + r 2 {s) so that we have (we omit the dependence on s) 

(2.29) - (rf + r| - 2rir2 cos(0i - ^ 2 )) 

/ Q 

= 2 rir 2 (l + cos(0i - ^ 2 )) = 4rir2 cos^ ' ^ 


The claim thus follows. 

The wave kernel on ( 74 ,^ is therefore given by 0, 1 or 2 times this quantity, 
according as there are 0 , 1 or 2 values of s > 0 satisfying t = ri(s) + r 2 {s), as 
discussed above. This agrees with the expression (2.1a)-(2.1c) obtained by Cheeger- 
Taylor. 


3. The microlocal structure of the wave propagator on Cq 

We now analyze the structure of the Schwartz kernel E of the sine propagator 
on the cone Ca of generic cone angle a. First let us recall the definitions of the 
geometric and diffractive Lagrangians and their intersection: the geometric (or 
“main”) Lagrangian is 

(3.1a) = rl +r 2 — 2 rir 2 cos(0i — ^ 2 ) and |0i — 02 | ^ t} , 

the diffractive Lagrangian is 

(3.1b) =^iV*{t2 = (ri+r2)"}, 
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and their intersection is the singular set 
( 3 . 1 c) 

In particular, we note that pr(I]) = = (ri + r2)^ and — di = ±7r}. 

To do this, we use Friedlander’s representation of the sine wave kernel on 
the cone of angle a, which expresses, in effect, this wave kernel as the a-periodized 
sine wave kernel on the cone of angle oo. Because of this, the wave kernels on two 
different cones Ga^ and Ca^ are closely related. We use this fact, together with our 
complete understanding of the case a = 47r from Section 2 , to prove the following 
theorem for any cone. 


Theorem 3.1. The Schwartz kernel E of the sine propagator W(t) on the Eu¬ 
clidean cone Ca is an intersecting Lagrangian distribution of class 


{rir2y 


r 




3 . 1 . Priedlander’s construction of the wave propagator. To start our study 
of the sine propagator near the singular set S, we recall Friedlander’s construction 
of the Schwartz kernel of W(t) from [FriSl]. 

Let G(j/, z) be the Lj^^-function on given by 

(H{y-\-cos{z))H{TT-\z\), y<l 

(3.2) G(j/,z)"=IW 

I- < arctan 

Form its periodization with respect to the map 

9 {y, z) I— > {y,z + a) & K^, 
and denote the resulting function by Ga{y, z); concretely. 


arccosh(j/) 


+ arctan 


arccosh(y) 


y>l. 


Ga{y, ^) = X! ^ • ^)- 

feGZ 


We may thus view Ga{y, z) as a function on K x (M/aZ). Now, define the operator 
A : K X (K/oZ) — > K X G° X G° as the composite A = A3A2A1, where 

• Al = [dy\ ^ is half-derivation in the y-variable, that is, the composition of 
differentiation in y with the fractional integral operator with kernel given 

by I2/-2/T®; 

• A2 = F* is pullback by the map 


F(t,ri,6»i;r2,6>2) 



• and A3 is multiplication by the factor 2Tr-^2rir 


Proposition 3.2 ([FriSl, Hil 05 ]). The operator A is a Fourier integral operator 
associated to the Lagrangian relation 

( 3 . 3 ) Af iV*|y = ~ ~ and z = 61- 62} , 

{ 2rir2 J 

and the Friedlander distribution AG^ on Kj x G° x G° is well-defined and equal to 
the Schwartz kernel of the sine propagator, E. 
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The important feature of Friedlander’s construction for us is the ease with which 
it decomposes E into pieces which are either associated to the geometric wave, the 
diffracted wave, or their intersection E. We use this to show the structure of E 
near S is the same (up to a purely diffractive term) for all cone angles a. 


Proposition 3.3. Let Cai and Ca^ be two Euclidean cones. There are isometric 
neighborhoods C x C°^ x C°^ and C x C°^ x C°^ of the set 

= (ri + r2)^ and 6i - 62 = ±7r} = pr(E) 

on which 

( 3 . 4 ) e/- 5 (F/,A'^), 

where Ea is the sine propagator kernel on Mt x C° x ( 7 °. The key point is that 
( 3 . 4 ) is purely diffractive. 


Proof. Let us start with the case where oi and 02 are both greater than 27 r. From 
Proposition 3.2 we know that Ea = AGa, so we may prove this proposition by 
showing an analogous statement for the periodized function Gq. We note that the 
projection of E to the base x G° x G° corresponds to 


y = 


f -rl- r| 
2rir2 


= 1 and z = 9 i — 62 = ±7r 


in the original (y, 2;)-coordinates used to define G. 

Let a* min(ai,a2)j and set e | (a* — 27 r). Choose a smooth bump func¬ 
tion p S C(?°(M2) satisfying p{z) = 1 when l^l < | and p = 0 when \z\ > e. From p 
and G we define the following: 


p ±( z ) p{zttt) 

p°{z) (1 - p+{z) - p-{z)) ■ l{|„|^^}(z) 

p^{z) (1 - p+(z) - p-{z)) ■ 1{|„|^^}(Z) 


G±(y,z)''=l'G(y,z)p±(z) 

G°(y,z) 1 l'G(y,z)p°(z) 

G“(y,z)"=l'G(y,z)p“(z). 


Thus, G = G° -I- G°° -I- G+ -I- G~. Using ( 3 . 3 ), and the calculus of wavefront sets, 
we see that these pieces of G correspond to the geometric wavefront, the diffracted 
wavefront, and a small neighbourhood of E, respectively, after periodization and 
the application of A. 

Now, consider the a^-periodizations G*^. (y, z) = of these 

distributions for j = I and 2 . By choosing e as above, so that aj > 27 r -|- 2 e, we 
have on the set {z S (—tt — e, tt -|- e) } 

Ga, (y, ^) - G„, (y, z) = G- (y, z) - G“ (y, z) 

since G^^ = G^^ ^ad G° ^ = G®^ here. Therefore, if we view Ea as the restriction 
of AGa to the fundamental domain for the periodization, we may set 

F~^{R X (±7r — £, ±7r -|- e)) and conclude 


E — E 

■*-^CXi ^0 l2 


v± 


= zl[G“ -G-] G J-J(U±,A°), 


since AG°° is purely diffractive. This establishes the result in the case oi, 02 > 27 r. 

Finally, to extend to general cone angles a we use the method of images: distri¬ 
butions on Ca may be represented as a-periodic distributions on its A^-fold cover 
Gtvoj where N is any positive integer. The result holds using Ejqa in place of Ea 
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by the above, and we may recover the result for Ea by restricting to a single period 
of length a in the angular variables. □ 


Proof of Theorem 3.1. Theorem 3.1 follows 
Proposition 3.3 . 

We conclude the microlocal structure of 
as a corollary of this result. 

Corollary 3.4. The Schwartz kernel U of 
K X (7° X (7° is an intersecting Lagrangian 

rfK~h-^R X (7° 

where (A°, A^) is the forward/backward part of the intersecting pair (A*^, A*^), i.e., 
the pair given by intersecting (A°, A^) with {(t,T) x T*C^ x T*Cf^ : =Ft > 0}. 

Proof. We know from Theorem 3.1 that W, the sine kernel, is in the class 

(7° X (7°; A'^, A^). By taking a derivative in t, we find that cost\/A is in the class 

/-i/4(]^ X Cq X (7a; A*^, A*^). We can write 

cos t^/A = ^ . 

Since is annihilated by the operator {Dt ± v^), which has symbol r ± |^|, 

we see that its wavefront set is contained in {=Fr > 0}. Therefore, is 

microlocally identical to 2 cos t'/A on A°, and microlocally trivial on A^. □ 

Remark 3.5. In [Hil05], a similar argument is used to pass from the sine kernel 
to the half-wave kernel but the factor of 2 has been incorrectly omitted. 


immediately from Proposition 2.8 and 

□ 


the half-wave kernel U K, 




the half-wave group U{t) = e 
distribution in the class 


xC':;A°,A^), 


Example 3.6. Starting from expression (2.26), this procedure yields the following 
expression. On the cone of angle 47 r, for 9i close to tt and 02 close to 0 we have: 


Ui^{t,qi,q2)m - ^ 


^OO poo 


0 ^0 


^i4){t,qi,q2,s,uj) 


Sin 


< ^l(s)+^ 2 (^) 


) 


(ri(s)r 2 (s))' 


Lv dsdijj\dqidq2 ^ , 


(where ~ means equal modulo (7°°) in which 0, rj(s), 9j{s) are defined as in (2.26). 


Remark 3.7. We emphasize that the novelty in Theorem 3.1 is the precise determi¬ 
nation of the structure of the wave kernel near the singular set E, the intersection 
between A*^ and A'^. Indeed, the Lagrangian structure of the wave kernel near 
A^ \ A"^ (where the cone point plays no role, due to finite speed of propagation) 
follows from classical work of Hormander [H 6 r 68 ] (also together with Duistermaat 
[DH72]). On the other hand, on metric cones (of any dimension), Cheeger and 
Taylor showed that the wavefront set of the wave kernel is contained in A*^ U A° 
and showed the Lagrangian structure of the wave kernel near A"^ \ A^ [CT82a, Sec¬ 
tion 2], [CT82b, Section 5]. More generally, on spaces with cone-like singularities, 
Melrose and Wunsch [MW04] proved that the wavefront set of the wave kernel is 
contained in A*^ U A°; morover, they also showed that the diffractive singularity 
is {n — 1)/2-order more regular than the geometric singularity. Notice that this 
difference in order agrees with our results, since for an intersecting Lagrangian dis¬ 
tribution, the order on Aq (here, the diffractive Lagrangian) is always smaller than 
the order on Ai (here, the geometric Lagrangian) by i. This also shows that our 
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result is restricted to dimension 2 : in higher dimensions, it cannot be true that 
the wave kernel on a cone is in the Melrose-Uhlmann calculus. In the latter case it 
would be interesting to know whether the kernel lies in the class of distributions that 
are constructed in [GU81] and that generalize the Melrose-Uhlmann construction. 

3.2. Proof of Theorem 0.3. Let us first consider the case of a diffractive geodesic 
of length t* joining to q* with a diffraction angle of tt. 

Remark 3.8. It may seem peculiar to use q 2 as the starting point and qi as the 
final point of the geodesic, but this is coherent with searching for an expression for 
U{t,qi,q2). 


We can use a Euclidean system of coordinates such that 

• ^2 corresponds to (—r 2 , 0 ), 

• the geodesic corresponds to the horizontal line starting from (—r 2 , 0 ). 

This Euclidean coordinate system can be uniquely extended to a local isometry^ 
from \ {(0, y), y > 0} into Cq. We will freely use this local isometry to identify 
points (qi, < 72 ) in a neighbourhood of ( 91 , 92 ) ''^ith their preimages in 

The point 9 * corresponds to (r*, 0) in this system of Euclidean coordinates. The 
geodesic between 92 and 91 is horizontal and it can be seen that it is geometrically 
diffractive with angle -|-7r since it is the limit of horizontal geodesics approaching 
from below. For s > 0, we denote by p+(s) the point with coordinates (0, —s) in 
this Euclidean system and we set 


(^+(t,9i,92,s,w) 


def 


92 -p+(s)| -f |9i 


p+{s) \ -t 


w. 


where 1 9 — 9 ' | denotes the Euclidean distance in . 

When the angle of diffraction is —tt we can proceed similarly. The diffractive 
geodesic is now the limit of horizontal geodesics from above and the cut is now 
{(0,y), y < 0}. We then define P-{s) (0,s) and 


91 , 92 ) 


I92 -P-(s)| + | 9 i -P-(s)| -t 


' (jj 


Lemma 3 . 9 . In either situation, locally near (U, 9^,92) G K x ( 7 ° x ( 7 °, (j)± is a 
phase function for the intersecting pair (A^, A^). 

According to corollary 3.4 and to section 1 there exists a symbol Oq, such that, 
locally near (U, 9 i, 92 ), we have the expression 

t^a(L9i,92) = (27r)“^ / / e*‘^±(*’'^i’'*"’®’‘^)aQ,±(i,9i,92,s,a;)dsdw|d9id92|U 

s>0 J u;>0 

Moreover, aa,± has an asymptotic expansion of the form 
(3.5) tta 92 , 5 ) 0 ;^“^. 

fc >0 

The only thing left to prove is the relation with the geometric theory of diffrac¬ 
tion. This is done by computing the leading amplitude of Ua near the diffracted 
front and away from S, and comparing it with Proposition B.l in the Appendix. 


4 


If a > 27 r this isometry is actually one-to-one onto its range. 
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Starting from the preceding expression and using the methods and results of 
section 1 , the leading term on the diffracted front is given by 


Ua{t,qi,q2) -(27r) ^ f 

J a;>0 

We compute ds4>±{t,qi,q2,s 
equation B.13. We obtain 


= 0 , 0 ;) = ±(sin 0 i +sin 02 )a;, and compare 


«a,±,i(gDg2,5 = 0) 
±i( sin0i + sin02) 


27r(rir2) ^ S'a(6»i - 6»2) 


so that finally 

(3.6) aa,±,i(t, 91 , 92 , s = 0,o;) ~ ■ 

This is the last statement in Theorem 0.3. 


Sc.{ei-92) 

(rir2)5 


• [sin 01 + sin 02 ] • w. 


Remark 3.10. This formula actually gives a way of computing Sa if we know the 
symbol in the Melrose-Uhlmann representation. For instance, starting from the 
formula in example 3.6 for the propagator near a diffractive geodesic with an angle 
—TT on a cone of angle Itt we derive 

_ 1 

a 477 ,-,i( 9 i, 92, s) = -i (ri(s)r 2 (s)) ^ sin 
The preceding formula thus yields 


^ 0i(s) + 02(s) ^ 


SiTr{di — 02 ) 


-27r 

—i (sin 01 + sin 02 ) 



— 1 sin (^(01 + 02 )) 
27r (sin 01 + sin 02 ) 



This agrees with the formula (B.12) in Appendix B. 


^ 0 l(O)+ 02 (O) ^ 


Remark 3.11. It is interesting to note that aQ,±,i( 9 i, 92 , s = 0,a;) is actually a 
regularization of the symbol on the diffracted front. The latter blows up when 
approaching the intersection and this formula gives an effective way of regularizing 
the contribution of a diffractive geodesic when the diffraction angle approaches ±7r 
(compare with the approach of [BPSOO]). 


4. The wave kernel after two geometric diffractions 


Theorem 0.3 can be used to understand the half-wave propagator on a ESCS 
after microlocalization along a particular diffractive geodesic. We will now present 
two applications of this method. A systematic study leading to a better knowledge 
of wave-invariants of a ESCS will be done elsewhere. 

Now that we have the basic structure of the half-wave kernel on the cone Cq,, 
we next determine the structure of the kernel after two diffractions on a Euclidean 
surface with conic singularities (ESCS). While one could continue to calculate the 
structure for an arbitrary number of diffractions and any kind of diffraction, we 
will focus on two geometric diffractions since this is the first case for which our 
approach yields a significant improvment on the existing literature. 
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Let X be an ESCS as described in the Section 0, and let q* and qj be two points 
in X with a geodesic 7 of length t* > 0 between them. Denote also by the 
covector in T*,X of the bicharacteristics that projects onto 7. 

Our aim is to find an oscillatory integral representation of the Schwartz kernel 
of the operator 


AiU{t) ^2, 

where Ai G is microlocalizing near (q*,^*) and U{t) g-*iv^ jg 

wave kernel at time t with t close to t*. 

In order to fix notations we assume the following. The geodesic starts at q^ then 
hits a cone point p2 then a cone point pi and finally ends at q^. We denote by a 
the distance (along this geodesic) from q^ to by 6 the distance between p2 and 
Pi and hy c = t* — {a+ b) the distance from pi to q^. Moreover, we suppose that 
this geodesic passes geometrically through both two cone points P2 and pi; i.e., 7 
is locally a limit of non-diffractive geodesics. 

Every geodesic with only one diffraction, which is geometric, is a limit of non- 
diffractive geodesics. For a general diffractive geodesic with several geometric 
diffractions, it may happen that, locally, the geodesic is a limit of non-diffractive 
geodesics, but not globally. However, in our case, since 7 has only two geomet¬ 
ric diffractions, it is always such a limit of non-diffractive geodesics (see [Hil06]). 
We show this by generalizing the construction we did for the geometric diffractive 
geodesic on a cone^. 

We start with a Euclidean coordinate system at q^ such that q^ corresponds to 
(—0,0) and the geodesic is horizontal and we try to extend this coordinate system. 
In the extended system the geodesic will correspond to the horizontal segment that 
joins (—0,0) to (6 -I- c, 0) so that p2 will correspond to (0,0) and pi to (6,0). We 
remove from the cuts 


cut2 {(0, e2s), s > 0 }, 
cuti {(6, eis), s > 0 }, 


in which e^, i = 1,2 is such that the angle of diffraction at Pi is Ciir. Exploiting the 
flatness of X, the original coordinate system can be extended to a local isometry 
from an open set E C \ (cuti U cut2) that contains the horizontal segment. If 
both €i have the same sign then 7 is the limit of non-diffractive geodesics that pass 
above the two cone points (or below the two cone points). If the have opposite 
signs, 7 is a limit of non-diffractive geodesics that cross the horizontal line between 
the two cone points. This case is illustrated in Figure 4.1. 

Using this local isometry we can define, for (<71,52) near (qj‘,52)) tbe functions 
\qi — Pi\ and 151 — 521 to be the Euclidean distance in of the corresponding 
preimages. 


^This construction is the same as the rectangles with slits that are used in [Hil06] 
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CUt2 


92 


'P2 

- *■ - ; 

'■■■■ TT ■■■' 


7 





cuti 


Figure 4 . 1 . The geodesic 7 passing through two cone points, with 
a geometric diffraction in both cases 

For t close to t* we can then define the following Lagrangian submanifold in 
T*{X° X X°) 

A° N*{\q2-p2\+b+\pi-qi\=t], 

iV*{|g2 -Pil + \Pi - qi\ =t}, 

=*' A^*{|q2 -P2I + \P2 - qi\ = i }, 

A^ N*{\q2-q^\=t}^ 

It is straightforward to check that A^ corresponds to direct propagation, A^ corre¬ 
sponds to one diffraction at pi , A^ to one diffraction at p2 and A° to two diffractions 
in a row at P2 and pi. 

The aim of this section is the following 

Proposition 4.1. For t > 0 fixed near t*, the Schwartz kernel of AiU{t) A2 is an 
intersecting Lagrangian distribution of order 0 associated to the four Lagrangian 
submanifolds (Aq, Ai, A2, A3). 

Proof. We begin with a decomposition of Ai Li (t) A2 in which only one conic point 
plays a role in each factor. This is straightforward: we choose a time to G (oj a + b), 
say to = a + and we write AiU{t) A2 = {AiUft — ttfi)) {U{tQ) A2). In terms of 
their Schwartz kernels, this is 

( 4 . 1 ) tC[Ait^(t) A2](91,92) = [ K:[AiU{t-to)]{qi,q) ■IC[U{to)A2]{q,q2)dq. 

J X 

Due to the assumptions on the microlocalizers the only points 9 that contribute 
to the singularities of ( 4 . 1 ) are points near (&/ 2 , 0 ) in our coordinate system. In 
each factor of the composition above, the singularities of the half-wave kernel only 
meet one cone point. Thus, modulo smooth errors, we may replace the half-wave 
kernel by the half-wave kernel on an exact cone in each factor, allowing us to use 
the results of Section 2 . 
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More precisely, in ( 4 . 1 ), to obtain a singularity q2,^2) in the canonical re¬ 

lation of AiU{t) A2, we must have (q,<72 j ^2) in the canonical relation of W(to) ^2 
and (qi, ^1; q, i) in the canonical relation of AiU{t — tQ). For t sufficiently close to 
t*, this implies that q is close to the point (&/ 2 , 0 ). That is, up to a error, we 
may insert a cutoff function x^(q) into ( 4 . 1 ), where y is supported close to ( 0 , 6 / 2 ): 

( 4 . 2 ) ~ 

JC[AiU[t) A2]{qi,q2) = [ IC[AiU{t - to)] {qi,q) ■ x^{q) ' /C[6f(to) ^i](g, 92) dg 

Jx 

modulo errors. Moreover, restricting the microlocal supports of Ai and A2 if 
needed, we may assume that the support of x is contained in a ball that is isometric 
to the corresponding ball in By the above assumptions on the geodesic 7, the 
half-wave operators U{to) and U{t — to) in the compositions x(g)6f(to)^2 and 
AiU{t — to)x{q) can be replaced (up to a smooth error) by the corresponding 
wave kernels on the exact cones with cone points pi, resp. P2, which we know from 
Section 2 are intersecting Lagrangian distributions associated to the diffractive and 
main fronts. That is, we can express the Schwartz kernel of x(g)6/(to)^2 in the 
oscillatory integral form 

p(X> pOO 

( 4 . 3 ) ( 27 r)"^ / / e"‘^^^'>'‘>^'*°’''^’‘^'^'>a2{q,q2,to,S2,uJ2)duj2ds2 

Jo Jo 

where (j)2 is the phase function 

( 4 . 4 ) (j)2 [|g - P2(S2)| + |P2(S2) - g2| “ ^o] ’ ^2, 

where ^2(^2) has coordinates (0,-6252) and | • | denotes the Euclidean distance in 

Similarly, the Schwartz kernel of AiU{t — to)x(g) has the oscillatory integral 
representation 

pOO pOO 

( 4 . 5 ) ( 27 r )"2 / / e'‘'^^^‘>^'‘‘’*~*°’''^'‘^^'>ai{qi,q,t-to,Si,uJi)duJidsi 

Jo Jo 

where (/>i is the phase function 

( 4 . 6 ) (/i [|gi - pi(si)| -k |pi(si) - g| - (i - to)] ■ wi, 

where now pi(si) has coordinates (6, —eiSi). Here, Ui is smooth, supported in 
Wi ^ 1, and is a symbol of order 1 in Wi. 

Therefore, ( 4 . 2 ) is given by an oscillatory integral (up to smooth errors) of the 
form 

p p poo poo 

( 4 . 7 ) ( 27 r)"‘^ / / / / ai(gi,g,t - to,si,a;i) 

J X.q J Si—0 J S2 —0 

X 02(g, g2, ^0, 52, W2) dsi ds2 duJi du}2 dq 

• • d0f . . 

We now show that in the overall phase function $ = (j)i + 4>2 we can eliminate 

the variables (9,012). This is possible if the following non-degeneracy condition is 
satisfied: 

( 4 - 8 ) d(q,^,)$ = 0 ^ detd^q_^^) (q,^^)$^0. 

The condition dg^ = 0 implies that q is on the segment [p2(52),Pi(51)] and that 
oil = UJ2- The condition d^^^^ = 0 implies that q is at distance to — 1^2(52) — g2| 
from ^2(52)- Since the non-degeneracy condition is coordinate free and has to be 
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verified with fixed si, S2,uji,qi, q2 we can choose for q cartesian coordinates (x, y) 
in a rotated and translated coordinate frame, such that the origin corresponds 
to the critical point and the conical points Pi{si) have the following coordinates: 
Pi(si) = ( 5 , 0 ), P2{s2) = (—^jO) with positive A, B. We observe that A and B 
depend on all remaining variables. 

In these coordinates we have (we only keep (x, y, ^2) as variables since the other 
ones are fixed) 


<^{x,y,u}2) = wi • \qi -pi(si)| + \/(5 - a;)2 + _ [t - to) 


+ W2 


-/(X+^F+l/^ + _ p2(s2)| - to 


We compute that 

W * - (^ - (^ + 

" V(5 - x)2 + y2 7(A + cr)2+/’ 

( 4 . 9 ) 

^{B - x)^ + j/2 ^(^ A + x)'^ +2/^’ 

= v/(A + a ;)2 + y 2 4 . -^ 2 ( 52 )! - to 


The critical point is easily seen to be (x = 0 , y = 0 ,022 = wi). We can then compute 
the Hessian of <i> in the (cc, y, W2)“Variables and evaluate it at the critical point: 


( 4 . 10 ) 


0 1 
0 ujiC 0 

1 0 0 


„ def 1 1 

^ - A^B- 


The determinant is —Cuji < 0 so that the non-degeneracy condition is satisfied. 
It is straightforward to check that this matrix has two positive eigenvalues and one 
negative eigenvalue. The signature thus is 1 . 

Hence, using the argument of Hormander [H 6 r 71 ], we can write the oscillatory 
integral where we replace {x,y^uj2) by their values at the stationary point that we 
denote by qc- We obtain the oscillatory integral (writing w for wi) 

/OO pOO pOO 

/ / 51 ,^ 2 , Si, S 2 ,w)dsids 2 (ia;, 

-00 J si =0 J s 2 —^ 

where the phase function ^'(t, <71,92, Si, ^2, Wi) is seen to be 


(4 12) ^ [192-152(52)1 + |P2(S2) “ 9=1 + |9c -pi(si)| + |pi(si) - 9i| -t]wi 

= [192 -P 2 (S 2 )| + |P 2 (S 2 ) -Pl(si)| -b |75l(si) “ 9l I “ t]uJi, 

and the amplitude is given by 
( 4 . 13 ) 

a(t, 9i, 92, si, 52, w) = e*’"/‘^(a;C') ^/^ai(9i, 9c, t - to, 5i, w)a2(9c, 92, to, S2, w). 


We can now verify easily, using Definition 1.11 and Remark 1 . 13 , that 4 / pa¬ 
rametrizes the given system of four Lagrangian submanifolds. Indeed, a sim¬ 
ple computation shows that at (t*, 9)“, 92, w*) = (c, (c — a, 0), (—a, 0), 1) we have 
dt.j,si,s2'f'(t, 9 i, 92 ,‘^*, 0 , 0 ) = 0 . Moreover, explicit computation shows that at 
this point the differential d{d'^/dsi) is a nonzero multiple of dj/i, the differential 
d{d'^/ds2) is a nonzero multiple of dy2, and d{d'i'/doJi) has a nonzero dt compo¬ 
nent. Thus these differentials are linearly independent, implying that the localized 
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propagator is an intersecting Lagrangian distribution associated to the above sys¬ 
tem. It is not hard to check that the four Lagrangians correspond to no diffractions 
(A 3 ), one diffraction (Ai, A 2 ), arising from interaction with pi or p 2 respectively, 
and two diffractions (Aq). Finally, as a in (4.13) is a symbol in w of order 3/2, we 
see directly from Proposition 1.12 that the order of the distribution (that is, the 
order on A 3 , the direct front) is 0 (where t is treated as a parameter). □ 

To conclude this section, we compute the principal symbol of the wave kernel 
U{t) at the twice-diffracted Lagrangian Aq using Proposition 1.14. This 

amounts to computing d, to leading order in w, at si = S 2 = 0. We will do the 
computation in the case ei = —1 and 62 = -1-1, as in Figure 4.1. The other cases 
are similar. 

Clearly, from (4.13), we need the leading order behaviour of ai at Si = 0. This 
is given by (3.6). Substituting into (4.13), we find that when si = S 2 = 0, 

(4.14) 

d(t, 91,92,0, 0,w) = e*’"/'‘C'“^/^(27r)^5'„i(-7r - 6*1) sin did 91 ~ P^\ \Qc - Pi\) 

X S'a2(d2)sind2(|92 - P 2 I |9c -P 2 I) mod 

When Si = S 2 = 0 we have 

^ ^ A + B ^ \qc-pi\ + |9e -P 2 I ^ b 

AB |9c-pil |9c-P 2 I \qc - pi\\qc - P2\' 

Using coordinates where (ri,di) are polar coordinates for 9 ^ centered at Pi, we can 
simplify (4.14) to 

(4.15) d(t, 91,92,0, 0,w) = e*’"/'‘(27r)^S'Q,i(-7r-di)5'Q,2(6'2)sindisind2 

X mod ^1/2. 

Using Proposition 1.14, and the identities 

ff'si = = a;sin(di) and = wsin(d 2 ) 

\y\ |a^l 

valid when si = S 2 = 0 , we find that the principal symbol at the twice-diffracted 
Lagrangian Aq is 

(4.16) 


1 

'a{x,y,to,t,0,0,ujy 


i9(ri,di,d2,w,ri + b + r 2 -t) 

2 1 

27T 

-1 

(M 

Co 

d{x,y,uj) 



'TT 6 1 

1,1 Sa^ (.62) Sai (-TT - di) |driddidd 2 da;| ^. 
0 J2b2 


5. Contributions to the wave trace of an isolated orbit with two 

GEOMETRIC DIFFRACTIONS 

As a byproduct of our approach we can compute in a rather straightforward 
way the leading contribution to the wave trace of any kind of periodic orbit, thus 
generalizing [Hil 05 ]. We present here the case of an isolated periodic geodesic that 
has two geometric diffractions (and no other diffractions). 
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More precisely, we assume that the orbit 7 diffracts at pi and P 2 (not necessarily 
distinct) and that the angles of diffraction are —tt and + 7 r. We construct as before 
the rectangle with cuts that is associated with this periodic geodesic. We see that 
near pi the geodesic is locally the limit of non-diffractives geodesics that pass above 
Pi- Near p 2 it is locally the limit of non-diffractive geodesics that pass below. It 
follows that one cannot translate the orbit to a nearby periodic orbit, so that the 
orbit is isolated as a periodic orbit. 


Remark 5.1. If instead of translating the orbit we rotate it then we do obtain 
non-diffractive geodesics that converge to 7 on any interval [ 0 , T] but these won’t 
be periodic. 


Let q be a point on 7, we intend to compute 


(Tp{t) = TT[AlU{t)A2p), 


for t close to the period L, Ai is a microlocal projector near and p a bump 

function near q such that on the support of p the principal symbols of Ai and A 2 
are identically I on the lift of the geodesic. More precisely, we first choose Ai and 
A 2 such that for t close to L, any geodesic of length t whose starting point is in the 
microsupport of A 2 and whose endpoint is in the microsupport of Ai stays close to 
7 . The bump function is chosen afterwards. 

We construct the Euclidean system of coordinates as before: the periodic orbit 
lies along the x-axis, with cone points P 2 located at ( 0 , 0 ) and pi at ( 6 , 0 ), and we 
identify {x,y) with (x -I- L,y), where L is the period. 

According to Section 4, the Schwartz kernel of the half-wave operator AiU{t)A 2 
after two diffractions has the following oscillatory integral representation 

poo poo pOO 

(5.1) (27r)"®/2 / / / 

Jo Jo Jo 
where '0 is given by 

i^{t,qi,q 2 ,Sl,S 2 ,Uj) = [\q2-P2is2)\ + |P2(S2) -Pl(si)| -f |Pl(si) - qi\ - t] ■ U} 


and d is given by (4.15): 


a(t, 91,92,0,0,0;) ~ (27r)^ • 
We are thus lead to compute 
(5.2) ap{t) = {2Tr)-^/^ 


2 sin 6 *i • ^^(-Tr - 6 »i) • sin 02 • S'a 2 ( 6 ' 2 ) f 


OO pOO poo 


xjo Jo Jo 


(^16x2)3 


dsidS2dqduj, 


where we have set 92 = 9 and 91 = 9 -I- {L, 0). 

We choose to parametrize 9 by (x, y): the Euclidean coordinates near qj. In 
this oscillatory integral, we first perform a stationary phase with respect to y. We 
denote by yc the stationary (critical) point. We observe geometrically that (x,yc) 
is on the segment [pi(si),P 2 (s 2 )]- Moreover, we compute 

|5yV’(L (a;,yc),0,0,o;)| = |^|| 


It follows that the critical point remains non-degenerate for small (si,S2)- Since 
geometrically, it is obvious that the critical point is a minimum, it also follows that 
the signature is -|-1. 
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We now observe that the phase, when evaluated at the critical point becomes 
independent of the remaining x. More precisely it is given by ij} where we have set 


( 5 . 3 ) 




■v /62 + (si + S2)2 + \/{L- 6)2 + (si + S2)^ - t 


• CJ. 


We thus obtain after applying the stationary phase: 

pL poo pOO pOO 

/ / / e^^A{t,x, si, S2, a;)p((a;, j/c)) dsi ds2 duj dx, 

Jo Jo Jo Jo 

where A is a symbol that, at leading order and for Si = S2 = 0, reads 
A{t, X, 0 , 0 , w) - a(t, (x + L, y^), (x, y^), 0 , 0 , w) 


x \dy'tpit,{x,yc),0,0,uj)\ 2 p(x,?/c) 

^ ^ sin( 6 ii) ■ Sai (-TT -0i) ■ sin( 6 > 2 ) ■ Sg^ ( 6 * 2 ) ^ \/\x\\L - 6 - x| 
-y/fo • \x\\L — h — x\ y/L — b 

. sin0i • S'o,, (-TT - 0i) • sin02 ■ 5 'q2(^2) 

^ i ■ - , -• UJ. 

VKL - b) 

It remains to evaluate an oscillatory integral of the form 

pL poo poo poo 

/(t) / / / / X, si, S2, w) • w dsi (is2 da; 

Jo Jo Jo Jo 


UJ 


in which A is a symbol in w. 

If we forget the restriction on the domain for (si, S2), this is a standard oscillatory 
integral and the phase has a smooth submanifold of fixed point. The restriction on 
the domain makes it a little less standard. Although we could perform a general 
treatment for this kind of oscillatory integrals, in our case, the nature of the phase 
allows for a more direct computation. 

We first make the change of variables m = si + S2, x = si — S2. In these coor¬ 
dinates, ^ is independent of v; we write '0(t, u,a;) for the phase expressed in these 
coordinates. Notice that, by ( 5 . 3 ), it is a smooth function of u^, and is stationary 
in u only at u = 0 . The domain of integration becomes m > 0 and —u ^ v ^ u. We 
obtain the integral 


m = 



A{u, V, uj)dv 


du duj. 


Since the factor in square brackets vanishes at m = 0 , the leading contribution of 
this integral is obtained by performing an integration by parts in u. To do this we 
write 

•4j{t,U,Uj) = {t — L)uj + tp{t,u,uj), tp{t,u,uj) = Olu^), M —>■ 0 . 

Then we have lim„_>o = d'^'ijj{Q,uj) 7^ 0 . We obtain 

pOO 

I{t) ^ i / 0, w)(i 9 ^i^( 0 ,a;))“^dw 

Jo 


where the index 0 means that we have taken the principal part of A. 
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It remains to evaluate all the quantities in our case observing that when si,S2 
go to 0 , 9 i go to 0 and 02 go to tt. Using (B.ll) and ( 5 . 3 ), we have 

lim sin 01 (—tt — 0i) = — 

lim sm02Sa2(02) = -J- 
62^tt It: 

Putting everything together we obtain: 

ap{t) ~ • £ P(^. 0)dx 

^ 7 0)dx 

The contribution of the whole periodic orbit is obtained by using a covering 
argument (i.e. choosing carefully near each point ^1,^2 and p so that in the end 
is identically 1 in a neighbourhood of the geodesic). In the process, we have 
to be careful near the cone point. The contribution of a (small) neighbourhood of 
the cone point can be computed using the following trick (that is already used in 
[Hil 05 ] and [Wun 02 ] ). 

Suppose Pc is a function that is identically 1 near pi. We want to compute 

o’pc(0 Tr(Z^(t)pc) 

We insert microlocal cutoffs so that, up to a smooth remainder we have 
(Tpc(t) = TT{AiU{t-to)A2U{to)pc)- 
Using the cyclicity of the trace we need to calculate 

(Tpc(t) = TT{U(to)pcAiU{t - to)A2). 

In the latter expression, thanks to the cutoffs, all the operations (composition and 
taking the trace) take place away of the conical point. So we can proceed as before. 

In the end, if we sum all the contributions, it will amount to sum all the contri¬ 
butions / p{x, 0)dx and this will give the length of the geodesic. 

We obtain the following proposition. 


Proposition 5.2. On a ESCS, the leading contribution to the wave trace of an 
isolated periodic diffractive orbit with two geometric diffraction is 


1 


• a/5(L — b) ■ {t — L — iO) ^. 


Remark 5.3. As a point of comparison, we recall the analogous leading-order 
contribution of a nondegenerate closed orbit 7 on a compact, smooth manifold in 
the trace theorem of Duistermaat and Guillemin [DG 75 ]: 

{2TT)-^Li‘^-' \ld- P^\~^{t- L-iO)-^. 

Here, P-y is the Poincare return map in the directions transverse to the level set of 
the symbol and to the flow direction, and P'< is a Maslov factor (with the Morse 
index of the geodesic). The singularity we obtain here from an isolated periodic 
orbit with two geometric diffractions is thus of the same order. 
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We can also compare this with the singularity contributed by a non-geometric 
diffractive periodic orbit with one diffraction, as computed in [Hil 05 , Theorem 2 ]. 
This has leading singularity (t—L—and is hence one half order more regular. 
On the other hand, the singularity contributed by a cylinder of periodic geodesics is 
to leading order {t — L — , from op. cit. which is half an order more singular. 

Notice that a cylinder of periodic geodesics necessarily has geometrically diff¬ 
racted geodesics at its boundary. In the second article in this series, we intend to 
use the analysis of the present paper to compute higher order terms in the wave 
trace singularity arising from such a cylinder. 


Appendix A. Domains of operators and admissible asymptotics at the 

CONE POINT 


In the course of our construction of the wave propagators on (74^, we needed in¬ 
formation about the domains of operators related to the Laplace-Beltrami operator 
Ag. The first such result was a description of the domain of the adjoint operator 
A~ * 


Lemma A.l. Let p G C°^{{ 0 ,oo)r) be a smooth cutoff satisfying p = 1 for r ^ 1 
and p = 0 for r ^ 2 . Then the domain of (A^) as an unbounded operator on 
is 


(A.l) 


D* 


D © Spang 


p,plog(r),pr2 exp 


± 


i0 


,pr 2 


exp 



Proof. Using the symmetry of Ag, we may decompose 23 * as 


23* = 23 © Null(Ag — ax) © Null(Ag — 02 ) 


for any distinct ai and 02 lying outside the spectrum of Ag (cf. [RS 75 ]). Moreover, 
nonnegativity of Ag implies that it is sufficient to let aj = —/ 3 | for distinct choices 
of fdj. Thus, let us suppose that u is an element of Null(Ag + ( 3 “^), i.e., 

(A. 2 ) (Ag +/ 3 ^) u{r, 0 ) = - ^ [r^d'^ + rdr - (/ 3 ^r^ - d^)] u{r, 0 ) = 0 . 

By separating variables using the spectral projectors ( 2 . 14 ), we may rewrite m as a 
Fourier series of the form 


u{r, 0) 



^Ug(r)exp 

jGZ 



Then the quality (A. 2 ) implies the corresponding equality 


(A. 3 ) 


(Lg+/3^)zig(r)^Af 



r'^dl + rdr 




iiffr) = 0. 


Introducing the change of variables s = fdr into (A. 3 ), this differential equation 
becomes 


(A. 4 ) 



s^dl + sds 




UjiP ^s) = 0, 


o2 

which is the modified Bessel equation, up to the overall factor of —Thus, the 
Fourier coefficients Uj must be linear combinations of the modified bessel functions 


li (s) and (s). 
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Now, observe that the condition that our original function u is an element of 
L^{C4tt) forces each of the Fourier coefficients Uj{l 3 ~^s) to be elements of the func¬ 
tion space L^(( 0 , oo)s, sds). Indeed, the Fourier decomposition in 9 induces a fac¬ 
toring 

L^C4.)=f{Z;L^{{0,oo)r,rdr)), 

and our change of variables identifies L^(( 0 , oo)^, rdr) with L^(( 0 , oo)s, sds)- This 
implies that the only admissible solutions to (A. 4 ) are 

uoiP~^s) = Ko{s), U±4{j3~^s) = Ki{s), and = 0 for |j| ^ 2 . 

These are the only modified Bessel functions which are globally in T^(( 0 , oo)^, sds), 
as may be easily gleaned from their asymptotics as s —>■ 0 and s —)■ oo in [AS 64 ]. 
Hence, 


(A. 5 ) Null(Ag+/ 32 ) 


Spanc 


Ko{Pr),Ki{l 3 r)exp 


'iO' 

~2 


> All (^r) exp 



Let p S C°°(( 0 ,oo)r) be a cutoff as in the statement of the lemma, and observe 
that 

[I- p{r)\Ko{Pr) and [I - p{r)] K i{j 3 r) 

are both Schwartz in r and vanish at the cone point. This shows they are elements 
of 2), which in turn implies that 2)* is equal to 


D © Spanc pKo{l 3 ir),pKi (/ 3 ir)exp 




,pKo{l32r),pKi{P2r)exp 


4 


for any two distinct choices of Pj > 0 . Similarly, since 

1 

Aro(x) = — log^^^ — 7 + 0(a;) as a:: —>■ 0 and Ki{x) = e~^, 

where 7 is the Euler-Mascheroni constant and r(2;) is the F-function, we have that 


Spanc pKo{Pir), pKi{Pir) exp 


4 


,pKo{p2r),pKi{P2r) exp 





r, 

1 


exp 

±2 

, p r 2 exp 

±2 




> (mod S) . 


This concludes the proof. 


□ 


The other piece of information about domains we needed was the expansion of 
elements of the Friedrichs domain I?2 at the cone point given in Lemma 2 . 1 . We 
now prove this lemma. 


Proof of Lemma 2 . 1 . The Friedrichs domain I?2 is characterized as the subspace 
of 'D* which is included in the Dirichlet form domain associated to Ag, i.e., those 
distributions u which are bounded in 


Qa^ (u) = {u, m)c2 + {u, Agu) c2 • 

As the Dirichlet form domain is precisely iL^(C'4^), we may conclude from the 
description (A.l) of D* that 


T’2 = J) © Spanc < p,pr^ exp 


, p r 2 exp 


iO 

'4 























WAVE PROPAGATION ON AN ESCS. I 


41 


since these are the only elements of which are elements of The 

lemma follows. □ 

Appendix B. Geometric theory of diffraction 

In this appendix we proceed with a construction of the kernel of the wave propa¬ 
gator that allows to compute explicitly the symbol on both Lagrangians and A° 
away of their intersection. For the diffracted part, this is known in the literature 
as the geometric theory of diffraction [Kel58] and we provide an interpretation of 
this construction based on scattering of waves on the cone C^- 

At the direct front, the symbol is just as it is on Recall that, on the 
half-wave kernel as a distributional half-density is 

(B.l) 

Let ei be a unit vector in the plane pointing from x to y, and let ( 01 , 62 ) be an 
oriented orthonormal basis. We write ^ = wei -|- pe 2 . Then the integral can be 
written 

(B.2) dpduj\dxdy\ ’. 

Assume t > 0. Then there are stationary points on the line {p = 0,a; > 0}. We 
can integrate out p, and to leading order (that is, replacing the expression with the 
leading term in the stationary phase expansion at p = 0 ) we get 

(B.3) (27r)-i y'e*(l"-^l-‘)“x(w)e-T da;|dxdj/|5 

where x G C“(M) is zero for w < 1 and 1 for w > 2. Thus, the principal symbol of 
this distribution at N*{\x — y\ = t}, tor t fixed, is 

(B.4) e~^xi^) \dydsdu}\^ 

for s the arc length along the circle {\x — y\ =f} and w the cotangent variable dual 
to \x — y\— t. 

We now return to Ca, the cone of angle a, and we restrict our attention to t > 0. 
On the diffracted front A*^ and away from the direct front, the half-wave kernel 
U (t) takes the oscillatory integral form 

(B.5) J K{r,6;r\d';u})\rdrd9r'dr'd9'\^ . 

The amplitude K{r,9;r',9'-,uj) is a symbol of order 0 in w, as follows from the 
kernel U(t) being of order — | (for each fixed t) at A°. We consider how this part of 
the propagator acts on a particular initial condition. Consider the exact solution® 
to the half-wave equation given by 

(B. 6 ) (I J|H (Ar)x(A) d\) 

where x G C°°(M) is supported in [2, 00 ) and identically 1 for A > 4, say, and where 
ly = ^ for some integer i. This distribution (B.6) is conormal to {r = — t} for 

®This solution is an example of the “plane waves” arising out of Cheeger’s functional calculus. 
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t <0 and to {r = t} for t > 0, as can be seen by using the expansion for the Bessel 
function as its argument gets large: 


(B.7) 


J\u\{z) = 


COS 2 — 


i=o 


where aj G S are the homogeneous terms in the expansion in z with ao{z) = 

1. Therefore, up to a smooth error, the solution (B.6) has the form 


Av9 


for a G 


. The singularities for t < 0, say t = —t*, take the form 

and for t > 0, say t = +t*, 

(B.9) (27r)-!e*'^®y’ dX. 


(B.8) (27r)-5e' 


On the other hand, if we apply the wave kernel to the initial condi¬ 

tion (B.8) we obtain (B.9) up to smooth terms. The direct front, away from the 
diffracted wave, does not contribute for t > <*, and the singularities come purely 
from the diffracted front (except at the intersection, where 6 — O' = ±7r (mod a)). 
Away from the direct front, applying (B.5) to (B.8) gives us 






X r' dr' dO' doj 


and after applying stationary phase in the (r', a;)-variables we obtain 


(27r)-5 J e*(’'-‘*)^(Ar)-5y(A) |e*l'^l^/2g^V4y 


(rr') 5 X(r, 0; r',0'; uj)e"'^' dO' dX. 


This must yield (B.9) up to smooth terms. Therefore the leading-order part of the 
amplitude S, viewed as a Schwartz kernel in the {9, 0')-variables, maps e"'^ to the 
quantity 

(27r) 5 e-"/2g-*kK i . 

Hence, the principal part of K corresponds to the operator 

-i(27r)5(rr')“5e“*’"'^^. 


since is an eigenfunction of V A§i with eigenvalue \v\. 

We can compare this principal part to the absolute scattering matrix S{X) for 
the cone Cq. This is, by definition, the map from the “incoming boundary data” of 
generalized eigenfunctions of VAgi with eigenvalue A, to the “outgoing boundary 
data”. These are the coefficients of respectively in the expansions of 

the generalized eigenfunction as r —^ oo. By inspection of the generalized eigen¬ 
functions 


Ju{Xr)e"'<^, 
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and (B.7), we see that this operator is —ie . Hence, we obtain at leading 

order and provided 9 — O' ^ ±7r (mod a) 

K{r,9,r',9') ~ {2TT)^rr')-^3^(9 - 9') 

where Sa{9 — 9''f is the kernel of the absolute scattering matrix for the cone of 
angle a or, equivalently the kernel of —ie ™ 

The principal symbol of the diffracted wave is therefore the leading-order part 


of 


K{r,9-r',9'-uj)\drd9d9'duj\i 


d{r,9-,r',9';u;) 


d{x,y,uj) 


— iTTy/ As 


{9,9') ■ {rr')"^\drd9d9'duj\-^ (mod S'"^) , 


which after simplification is 

(B.IO) ^/^So,{9-9')\drd9d9'duj\-^. 


The distribution Sa can be computed using Fourier series. Indeed, since it is the 
kernel of — ie we have 

rv ( ^ 


kG'Zi 

—i 
a L 




^(TT-e) 


k>l 


E< 

fc>i 




(B.ll) 


—i 

a L 

—i 

a 

-1 


1 -h 

1-f 






1 — e“ 1 — e“ ^ - 


2isin(^(7r-6»)) 2isin (^(tt- h 6»)) 
sin( —) 

V O' 


2a sin (^(tt - 9)) sin (^(tt -h 9)) 
In the case a = 47r, this simplifies to 

1 -1 


(B.12) 


^4.(0) = ^ 


Stt sin(^)sin(^) 47r cos f' 

Summarizing this computation we have the following proposition. 

Proposition B.l. Microlocally near the diffracted front and away from S, 
the leading part of the half-wave kernel Ua on the cone of angle a is given by the 
following oscillatory integral (using polar coordinates) 


(B.13) Ua{t,qi,q2) 


J_ f „i{ri+r2-t) 


27r 


(rir2) ^Sa{9i- 92)dw\dqidq2\ 


'cj>0 


with 

(B.14) 


Sc.{0) = 


-1 


sin(2^^) 

\ a / 


2a sin (^(tt - 9)) sin (^(tt -t 9)) 

Remark B.2. This coincides with Theorem 4 in [HilOS] up to the factor 2 that as 


been omitted there. 


'We have used the invariance by rotation to write this kernel in this form. 
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